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Introduction 


During last decades quantum chromodynamics (QCD), as a part of the 
Standard Model, established itself as a working theory of the strong interac¬ 
tions. We believe that quarks and gluons are the basic objects that make np the 
hadronic matter. This theory has been very successful in predicting phenom¬ 
ena involving large momentum transfer. Unfortunately, its direct applications 
are limited to this particular case, when the coupling constant is small and 
perturbation theory is reliable. However, in the hadronic world, the coupling 
constant is of the order of one and the perturbation series do not converge. In 
this domain the discrete space-time formulation of quantum chromodynamics 
provides a non-perturbative tool for the calculation of the hadronic spectrum 
and basically every other observable from first principles. Lattice QCD was 
used to study the questions of confinement and chiral symmetry breaking. 

Lattice QCD is usually formulated on a square Euclidean space-time grid. 
No new parameters or field variables are involved in this formulation, hence, it 
retains the fundamental character of QCD. In fact, it can serve two purposes. 
First, it acts as a non-perturbative regularization scheme. At finite values of 
the lattice spacing a, which provides an ultra-violet cutoff there are no infinities 
(as long as lattice is finite itself). On the other hand it gives us an ability to do 
direct simulations of lattice QCD on the computer using methods analogous 
to those used in statistical mechanics. One can in principle calculate all cor¬ 
relation functions of hadronic operators and matrix elements of any operator 
between hadronic states. These simulations are based on a Monte Carlo inte¬ 
gration of the Euclidean path integral and have therefore both discretization 
and statistical errors. 

The starting point for lattice QCD is the partition function in Euclidean 
space-time 

Z = J VA^V^e~ s (1) 

where S is the QCD action 

5 = f d 4 x - i>M^j (2) 

Here A4 is Dirac operator, fermions are represented through Grassmann vari¬ 
ables 0 and 0; F IIU is the Yang-Mills held strength tensor for gauge fields A 
which are elements of the SU (3) algebra. Fermion fields can be then integrated 
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out exactly, leaving us with 


Z = / VAn det M.e 


-Sg 


(3) 


The partition function is now an integral over background gauge configurations. 
The fermionic contribution is contained only in the very non-local term det M.. 
The model defined so describes zero temperature physics. To calculate the 
thermodynamical quantities at finite temperature one has to make the time 
direction periodic and set its length to be the inverse temperature. We write 
then 

S = ^ dx 0 J d 3 x Qf^F'"' - (4) 

where Xo refers to the time coordinate and V is the physical volume of the 
system. The gauge fields have then periodic boundary conditions and fermion 
fields - anti-periodic ones. 

The fermionic determinant is the major complication in all simulations on 
the lattice, due to its non-locality and other technical details. The model 
without this term, so-called pure gauge theory, has been in the focus of the 
interest for many years. It was successfully simulated by several groups || 
and the equation of state and the phase diagram were successfully studied and 
extrapolated to continuum 

The full model, however, suffers from severe problems of both calculational 
and theoretical type. Despite significant progress during last years ||, the 
time estimated for the continuum extrapolation of the results from the direct 
simulations of any of the existing fermionic actions (aiming, e.g. at the deter¬ 
mination of the quark masses) varies from several to hundred teraflops-years. 
Thus one would want to have an indirect method of calculating the quantities 
of interest, which would require significantly less computer time. One of such 
methods, applicable at high temperatures, is the so-called dimensional reduc¬ 
tion (0-[n|). In general it is a method of turning (D-fT)-dimensional system 
at high temperature into the D-dimensional one. Using the fact that at high 
temperature the domain of integration in (0) shrinks one can perturbatively 
get rid of the non-static time-like degrees of freedom. Such an effective model 
is much cheaper to simulate and is supposed to reproduce the long-distance 
behavior of the model. This method was successfully applied to QCD at high 
temperature and to the study of the electro-weak transition. At the same time 
the domain of its validity is not well-studied, and this is the primary aim of 
the current work. 

We focus on the (2+l)-dimensional pure gauge field theory, which is then 
reduced to a two-dimensional one. The source model does not include the dy¬ 
namic fermions and can be simulated with a high precision using the present 
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computer power. While being different from the standard (3+l)-dimensional 
system it shares with it many important properties, like static quark confine¬ 
ment; the infrared singularities are even more severe. The reduced model can 
be simulated even on personal computers. 

Current work is organized as follows. In the first chapter we define our 
model and operators of interest and review simple models of the dimensional 
reduction. In the second one we present the dimensionally reduced model and 
compare it to the original one. The third chapter is dedicated to the reduced 
model per se: its spectrum, phase structure and other properties. The fourth 
chapter will address the restoration of the Z( 3) symmetry and is followed by 
conclusions and outlook. 
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Chapter 1 

Reducing the degrees of freedom 


1.1 The Model 


Pure gauge chromodynamics on the lattice may be formulated in terms of closed 
loops on the square lattice A with link size a. We will denote its extensions 
in space and in time directions by L s and L 0 correspondingly. The lattice is 
periodic in all directions and its length in the time direction equals the inverse 
temperature 

L 0 a = l/T (1.1) 

On each link of the lattice we put an SU( 3) matrix U(x ; +); it is linked to the 
gauge fields A /t (x) by 

U(x] /i) = exp iA^x) (1.2) 

where go is the bare strong coupling constant. The main requirement for the 
action is the invariance under local gauge transformation generated by SU (3) 
matrices G(x) 

U(x; fi) — > G\x)U(x] /i)G{x + +). (1.3) 

The simplest such action is the Wilson one-plaquette action 

S(U)=h Y. (l-gFterrf/nu) (1-4) 


Here 



(1.5) 


is the lattice coupling constant. In (2+1) dimensions has a dimension and 
sets the scale of the model. The coupling /3s is inversly proportional to the 
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lattice spacing and thus proportional to the temperature (for constant Lq). 
The plaquette variable 

Up U (x) = U(x; n)U{x + //; v)U{x + /j,)^ 1 U(x; v)~ l (1.6) 

represents a loop around one plaquette. It is easy to show that in so-called 
’’naive” continuum limit when we send the link length to zero and expand the 
U fields around unity one recovers the Yang-Mills action 

S = \tt I d A xF (1.7) 

The partition function is 

Z— f dUexp(-S), (1.8) 


where dU is the invariant Haar measure. The expectation value of any operator 
is defined thus as 


f dUO(U) exp (-S) 
f dU exp (-S) 


(1.9) 


The operators which one frequently is interested in, and which will be the 
major focus of our consideration are the Polyakov loop, which is a trace of the 
product of the link variables along the time axis 


Lo—l 

L (x) — -Tr U(x, ia; 6), (1.10) 

3 *=o 

and the spatial Wilson loop, which is a trace of the product of the spatial 
link variables along the closed path. We will consider only rectangular Wilson 
loops, so for a rectangle R = [R \, R 2 \ 

W(R u R 2 ) = ^l[U(x,l) (1.11) 

x,l^R 


In (2+1) dimension this model has a deconfinement phase transition of the 
second order. The order parameter (on infinite size lattices) is the Polyakov 
loop. The critical coupling f3 c has been determined precisely in ||T^[ and is 
14.74(5) for L 0 = 4. We will later use this value as the basic unit and express 
temperature in units of the critical temperature, using for fixed L 0 


T (3 (3 

T c (3 C 14.74' 


( 1 . 12 ) 


Above the critical point there exists a finite string tension and the Debye 
screening mass is not zero; screening mass vanishes at the critical temperature. 
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1.2 Naive reduction 

In this section we will briefly review so-called “naive” dimensional reduction. 
It is a rough approximation of the full reduced model, which we will dis¬ 
cuss in Chapter 2. In fact for very high temperatures one may integrate out 
time-related degrees of freedom. Then we are left with (in our case) a two- 
dimensional lattice A 2 with Wilson action 

S 2 w = -f3 2 J2l ReTrU P (!- 13 ) 

x€A 2 

with the coupling /? 2 being L 0 times bigger than the original coupling of (2+1 )D 
model. 


@2 — L 0 f3 3 


(1.14) 


This approximation is particularly interesting due to the fact that this 
model can be analytically solved. The solution comes up to Gross and Witten 
39| , here we will sketch the main ideas and present the actual calculation of 


the string tension (we will need this result in Chapter 2) in Appendix 2. 

Let us consider an infinite (or with open boundary conditions) 2D square 
lattice as described above. Due to the nature of our 2D model we will use 
symbolics 1 and 2 for directions x and y, leaving 0 for the time direction. For 
this model the following gauge may be introduced. First, we set an arbitrary 
gauge matrix in the coordinate system origin and gauge all link variables in 2 
direction to one. Then, starting from this line, we set all links in 1 direction 
to unity (like in the static gauge, although we do not have time-direction now) 
and have thus: 


U(x, i) = 1, 
f/(0i + u2,2) = 1 


Now we introduce the following change of variables 

U(x+ i, 2) = W(x)U(x,2), 

so that the action factorizes and the statistical sum becomes 

J+lReTrW, = 




(1.15) 


(1.16) 

(1.17) 
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z = 


(1.18) 


dW exp 


(3- Re TrW 
o 


being a one-link integral. 

With the same argumentation one can get for the spatial Wilson loop with 
extensions R and T 

W L = (w) RT (1.19) 

with 


w — I dW-TrW exp 


The string tension is then 


/ 3-ReTrW 
3 


d 


<y = - In w = - log —z 


( 1 . 20 ) 


( 1 . 21 ) 


One can immediately notice that w is always smaller than one; thus string 
tension is positive, and one is in the confined phase for all coupling values. 
This is the property which all existing dimensionally reduced models have - no 
deconfinement phase transition. This reduction is useful only as a reference 
point; it is not possible to define Polyakov loops and the screening masses in 
this model. 


1.3 Classical reduction 

The next step in our consideration is so-called classical reduction. The name 
’’classical” refers to it being performed on the level of action, not taking into 
account the quantum fluctuations. Consider the general (we will do this in 
continuum, and use different normalisation of the potential to have the coupling 
constant in front of the integral) action for a pure gauge SU( 3) Yang-Mills 
theory in (2+l)D: 


S 



2 

d 2 x ^ TrF^F^ 

ft,v=0 


with the held strength thensor being 


F^iu = d»A v - d v A tl - ig[A v , A^ 


( 1 . 22 ) 


(1.23) 


For very high temperature it seems completely legitimate to ignore the time- 
dependence of the fields 


d 0 A u = 0 


( 1 . 24 ) 



Now we can rewrite F^ v in the following symbolic way (latin-style indexes go 
over spatial directions only): 


Fkm F o -do] T d m A() ig[A m , dp] 

— Fkm F -Dfcdo + -D m do 

with the covariant derivative being 

-Dfcd 0 = <9fcd 0 — ig[Ak, d 0 ] 

The action then takes the form 


(1.25) 


i r 

TrF km F krn + 2 ^ (D 

with the new coupling 


fcd 0 )" 


k,m= 1 


k,m= 1 


G 2 = g 2 T 


(1.26) 


(1.27) 


(1.28) 


This is the action with the same gauge group in two dimensions, coupled to 
the scalar held d 0 which transforms under the adjoint representation of the 
gauge group. All the fields are called static modes because they do not have 
any dependence on T. The latter is present only in the coupling constant. 

However, as we mentioned before, this procedure is classical and ignores 
the quantum fluctuations in the system. Detailed consideration performed in 
e.g. jni shows that the complete decoupling of the non-static modes is not 
supposed to be happening in this model. 
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Chapter 2 

Dimensional reduction 


2.1 Perturbation Theory 


Now we will turn to the actual derivation of the reduced model. We are inter¬ 
ested in gauge invariant operators, such as the Polyakov loop correlations and 
the spatial Wilson loops, and those may be computed in any gauge. We will 
therefore choose the following particular static gauge, which can be realized on 
the periodic lattice by first setting all A 0 to be independent of the imaginary 
time coordinate x$ 

A 0 (x 0 ,x) = A 0 (x) (2.1) 

This automatically makes the Polyakov loops static operators. Remnants 
of the gauge freedom are eliminated by adding the so-called Landau constraint 


2 

Y Y [A (x) - Ai(x- ai )] = 0 (2.2) 

xo ?=1 


Together with m this condition is called Static Time-Averaged Landau 
Gauge (STALG) 


D,[ 


Spatial Wilson loops are not static operators in 


this gauge. 

One can split the held in the static and the non-static components, Af atlc 
and A™ s which are defined by 


A i (x)=A? atic (x)+Ar(x 0 ,f) (2.3) 

XI ^rUo, 50 = o (2.4) 

XQ 


The non-static modes should be integrated out, leaving only the static compo¬ 
nents. We will do this in a perturbative fashion and will omit the superscript 
’’static” then. 
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It is well known that at the high temperature Z( 3) symmetry of the Wilson 
action is broken and the average of the Polyakov loop is not zero and is close to 
one of the cubic roots of one. One can eliminate this three-fold degeneracy by 
rotating the phase away and keeping all gauge variables close to one. In this 
case one can perform a perturbative expansion, using the Taylor series around 
unity in arbitrary parametrization satisfying ( p.2|) and ( |2.1| ). 

We obtained the effective two-dimensional action according to the standard 
perturbative techniques (1,0, [ |I0|). In momentum space the effective action 
consists of the classically reduced 3D action with all fields not depending on 
the 0 th coordinate (static modes only) and the one-point irreducible Feynman 
graphs with n— external legs and non-static internal lines, which we will denote 
by We have for the n — th order 


sfff (pi-Pn-l) = So n) (Pl-Pn-l) - r i n s } (Pl-Pn-l) ■ 




(2.5) 


It may be proven that f ^ are analytic in momenta due to internal lines having 
no singularities, so we can expand them around zero momentum and keep only 
the high-temperature contributions. 

Let us see which diagrams we need to calculate. We will limit ourselves to 
the one-loop approximation. For superrenormalizable theory in 3 dimensions 




( 2 . 6 ) 


where g r is the renormalized coupling, go - the bare coupling and A is the 
momentum cutoff parameter. Higher order contributions to R get higher 
powers of A, thus smaller powers of T. For diagrams each external leg brings 

1 /9 

in (g R A) ~ in the nominator. 

The diagrams with three external legs, which would lead to the terms pro¬ 
portional to the TtAq, are not permitted by the time reversal symmetry of the 
original action. This symmetry is discussed more in detail in chapters 3 and 4. 

The terms with four external legs are of the type f d 3 k/k 4 and thus UV- 
hnite. The restriction to the non-static modes k 0 ^ 0 makes them also IR-hnite. 
Therefore we computed the corresponding diagrams in the continuum, which is 
simpler and avoids cutoff effects. The Feynman rules were taken from ||6]] . The 


direct calculation was compared to the adaptation of the Nadkarni calculation 
0 to the D = 3 case. 

For the term proportional to the square of the ’’liiggs” held we need the 
H-(0), which is the one loop contribution of the non-static modes to vacuum 
polarization at zero momenta. Two Feynman graphs contribute, both linearly 
divergent with A. We computed them on the lattice (see Appendix 3), which 
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provides the cut-off of the order of 



a 


TL 


o 


(2.7) 


at high temperature. 


2.2 Action 

The resulting action depends thus on the parameters of the original theory L$ 
and /?3 (we will set a = 1 from now on). We make a change of variables 

A) (x) = ( 2 - 8 ) 

which will be useful later and express as the function of this field and the 
two-dimensional gauge field U. It consists of the following three parts: 

• pure gauge part S fy, which we had already in ’’naive” reduction 


= /3 2 V TV ( 1 - - Re U(x; l)U(x + ai; 2)U(x + a2; i)” 1 f/(a;; 2)" 1 ] 
xdh.2 ^ 3 ' 

(2.9) 

P 2 = L 0 P 3 ( 2 - 10 ) 

• gauge covariant kinetic term, depending both on U and 0, obtained by 
expanding the to the second order in A 0 (we had it already in ’’clas¬ 
sical” reduction) 


Tr (D t (U) p(x)) 

x £= 1,2 

with Di being covariant derivative 


( 2 . 11 ) 


Di (U) p(x) = U(x; i)p(x + ai)U(x ; t) — p(x 


( 2 . 12 ) 


Now one can easily see the reason for the normalization of the Higgs-held 
(|2.8|) - the classical limit of the kinetic term is Tr (dip) 2 . 
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• self-interaction term S^, which we will also call ” Higgs potential' 


S <f> = h 2 J2 Tr(/> 2 {x) + h 4 {Tr(j) 2 (x)) 2 

X X 


with the couplings 


h 2 

h 4 


9 

ttLo/33^' 

(l"K Lo + | log 2 - 1^ , 
9 

16ir/3j’ 


(2.13) 


(2.14) 

(2.15) 

(2.16) 


We note here that the second term in the action contains both terms of the 
type ( Tr(jr(x )) and Tr(jr{x) which are related to each other via 


(Trd> 2 (x)) 2 = 2 Tr<p‘(x) 


(2.17) 


in the case of 577(3) group (or, in our case, for the variables which are related 
to it via multiplication by a constant factor). 

The coupling h 4 comes from the diagrams with four external legs, while the 
quadratic coupling is expressed via 


ho = 


6 


nffito) 


(2.18) 


which can be numerically calculated from Eq.( A.52 ). We used, however, the 
scaling form of this quantity ( j2.14| ), that is its asymptotic value at large L 0 for 
the reason that we assume (and prove) that we work in the scaling region of 
the model (see subsection [2.4.2|) . The details of the derivation can be found in 
Appendix 3. 

Let us now have a closer look at this action. Llnlikc the quadric coupling 
/r 4 , which is positive and insures boundedness of the partition function at large 
fields, the quadratic coupling h 2 is negative. The potential thus has a shape 
typical for the gauge symmetry breaking by Higgs mechanism. Of course, this 
symmetry breaking is not expected to happen in two dimensions. In fact, the 
form of this coupling can be seen as the counter-term for the logarithmic ultra¬ 
violet divergence of the two-dimensional model, here prescribed in advance by 
the UV regularization of the original model. The actual infrared behavior of 
the model is highly non-perturbative. 


14 









In the static gauge the Polyakov loop operator is given in 2 and (2+1) di¬ 
mensions by the same function of the A 0 held. According to our normalization 
of the 0—held it reads 


L(x) = ^TV exp (i<j>{x) / \/t) (2-19) 

o 

with 

A = A = ! 

6 T 0 6 L 2 0 g\ 

being a dimensionless temperature. 

Along the line of constant physics (r = const ) the operator used to probe 
dimensional reduction remains unchanged, and it is thus consistent to compare 
the correlations of the Polyakov loops measured in three dimensions at the same 
hxed r value. The continuum limit is To oo, and large distance behavior 
can be checked within either one of the two models. 

Let us now consider the effective action S 2 ^ as a function of the dimension¬ 
less temperature r and the parameter of the original lattice Tq. After setting 



2 2 
a g 2 


6 

% 


( 2 . 21 ) 


and using the equations above we find the following form for the three dimen¬ 
sionless couplings 


2 2 
a g 2 


1 


tL\ 


( 2 . 22 ) 


_ 3 c 0 log T 0 + ci 

h2 “ 2vr tL\ 

/ J_LI 

4 647T tTq r 

One can easily see that apart from the logarithmic corrections all three cou¬ 
plings scale as 1 /Tq with To- As a function of r, h 4 is r-times smaller than the 
two other ones. 

Let us now write down the effective Lagrangian T e // as a result of the small 
a expansion of the effective action S^j. We make a substitution Ai —> , 

Y^ x —^ c+ 2 / d 2 x and use T and g 2 as the parameters instead of 03 and T 0 and 


(2.24) 
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take the limit a —> 0 everywhere but in the logarithmic term. We get 


c <fi = iE^+ Tr i- D « 


12 . 92 92 rji , 4 I r 

+-T^r-r^Trip + L C t 


c= 1 


32t r T 2 


Di<t> = di<j) + [>!*, </>] 

Fij = diAj — djAi + \Ai, Aj] 


r. _ 3^2 

CcT ~ 27 


2 r 


— log aT + - log 2 — 1 


Tr0 2 


(2.25) 

(2.26) 

(2.27) 

(2.28) 


This is a two-dimensional SU(3) gauge invariant Lagrangian for the adjoint 
scalar held 0; it is far from being the most general one. The gauge coupling 
#2 has canonical dimension one in energy and sets the scale. The non-kinetic 
quadratic term is the counterterm, suited to a lattice UV regularization with 
spacing a . The quadric self-interuction 

k = 9? f| (2.29) 

goes to zero as the temperature approaches infinity and presumably plays a 
marginal role. 

It is easy to identify the symmetries of the reduced action. The space 
inversion symmetry of the original model remains untouched, as well does 
the rotational symmetry in the space-like plane. The time direction reversal 
symmetry, as we pointed out above, turns now to the invariance of the action 
with respect to changing the sign of the higgs fields, so we define 


R t : (j) —> — <f> 


(2.30) 


2.3 Simulations 


To fulfill our aim of direct and detailed comparison, we had to simulate both the 
three dimensional and the two dimensional models. For the three dimensional 
part we used the Quadrics QH2 parallel computer of the Bielefeld University. 
We refer to [|28[ for the detailed description of this machine. The actual pre¬ 
sented data are combined data from C.Legeland and M.Luetgemeyer, which 
were crosschecked with the data obtained by the program independently writ¬ 
ten by the author (based on a different simulation design model of the Bielefeld 
Lattice Group). Both 3D programs use the Cabbibo-Marinari pseudo-heatbath 
update algorithm, based on updating the SU (2) subgroups of the 577(3) matrix 
using Kennedy-Pendelton algorithm. No additional technique like smearing 
was used. 
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The two-dimensional project was realized on our farm of PCs. It consists of 
one server and 10 diskless clients, each having AMD K6 processor at 300MHz 
clock speed. The cluster was gradually upgraded to a set of AMD Athlon 
processors, first to 600MHz and then to 1200MHz ones, which are (apart from 
the linear scaling with the clock speed) considerably more effective in floating 
point operations. Booting of the clients is done via traditional etherboot [[HJ 
which resides on the network cards and boots system from the server. The farm 
runs under Linux operating system from RedHat with appropriately modified 
kernel. It provides Message Passing Interface (MPI) services via Los-Alamos 
MPI Daemon (LAMD) from the University of Notre-Dame jp2 


The two-dimensional program (written in collaboration with P.Bialas using 
parts of the standard Bielefeld C++ Monte-Carlo code) consists of two main 
parts - master and slave. 


Master part: 


• runs on server 

• reads and parses the configuration hie 

• gives an appropriate set of parameters to all slave notes. Parameters can 
be either different for exploring a certain area in coupling constants or 
the same for accumulating statistics. 


Slave part: 

• runs on clients 

• gets the parameters 

• does the actual simulation 

• outputs the results 


The program is realized in an object-oriented way using classes lattice 
(SU(3) group) and higgs (SU(3) algebra). Measurements are done via classes 
basic_observables (all local observables), corr_observables (all correlations) and 
swLobseruables (non-local observables). All measurement classes are inherited 
from the class observable which does pre- and post-processing as well as output. 

For updating gauge fields we used the same algorithm as above, but after 
each subgroup update the resulting matrix was subjected to the Metropolis 
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question on the hopping term. This approach results in an acceptance ratio of 
about 95%. 

Alternatively we tried the multi-hit Metropolis algorithm. An update of 
a single link variable was attempted several times (typically eight) before at¬ 
tempting to update the next variable. The trial matrices were generated ac¬ 
cording to 

U — exp(iSA)U 

where 5A was a Hermitian matrix generated from the distribution 
exp(— eTrA 2 ). The parameter e was tuned to obtain the acceptance ratio close 
to 50%. Because generation of the trial matrices is expensive in CPU time 
(for it requires matrix exponentiation), we stored the values in the table for 
later use (that is we store both the matrix and its conjugate to ensure the de¬ 
tailed balance). This table was updated every five sweeps. The autocorrelation 
times were comparable to heat-bath but the CPU time required to perform the 
metropolis update was approximately two times bigger. 

The scalar fields were updated using multi-hit metropolis. At every point 
we tried as before typically eight updates. New scalar fields were obtained from 

A -> A + 8 A 


where 5A was generated the same way as above, using the same desired accep¬ 
tance rate. One sweep over lattice consisted of one update of the gauge fields 
followed by one update of the scalar fields. 

Operators we measured were: 


• the ” control” variables trA 2 , trA 3 and trA A 


• plaquette variables averages 


Polyakov loops (defined by (|1.10| ) in 3D and (|2.19| ) in 2D) and their 
on-axis correlations 


L(r) = ^ (L(x)L(x + r) + L{x + r)L[x)) — ( L(x )) 2 (2.31) 


where the averaging is performed over lattice and computer time. When 
necessary we will use indexes D = 2 + 1 or 2 specifying the dimension 
considered. 


• Spatial Wilson loops W(r\,r 2 ) (defined in the same way both in 2D and 
3D ) 


We used the following tests to test the program: 
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• Gauge Sector: remove interaction (kinetic) term in the update routine. 
Then the model is reduced to the pure gauge theory with action (|1.13|) 
and can be compared to the analytic solution. 

• Higgs Sector: set the f3 3 and the quadric coupling to zero. The model is 
reduced to the free-held model for the scalar fields which can be exactly 
solved. 

• Polyakov loop Correlations: same as in latter case, may be expanded 
and calculated in the small held limit 

The analytic solution and results for the two last test are presented in the 
Appendices 4 and 5. 

All runs in the main simulation of the two dimensional model were per¬ 
formed on 32 x 32 lattices with the parameter Lq set to 4. The dynamics of 
the gauge and Higgs sector turned out to be quite different. The integrated 
autocorrelation time for the plaquette was typically of the order of one while 
for the Higgs related local variables (e.g. trA 2 ) it varied between 100 and 400 
sweeps when /3 3 varied between 21.0 and 346.0. For each (3 value we collected 
typically one million sweeps. 


2.4 Polyakov loop correlations and Masses 


2.4.1 Comparison with the original model 


The Polyakov loops correlations defined as in Eq.( |2.31| ) were measured every 
ten sweeps. While due to the large autocorrelation times this was in a sense 
an ” oversampling” - it was not too costly in CPU time. The measurements 
were blocked by 50 and written out. Afterwards the data were analyzed using 
the double jack-knife method with typically 25 blocks. First, the jack-knife 
procedure was used in a standard way to get the values of the L{r) and the 
error estimates. Then these were used to obtain the fit (see below) parameters 
and the errors (which were uniform for all blocks of the second jack-knife) 
using the second jack-knife. For checking consistency we also tried a different 
number of blocks, and the errors, as expected, did not show any significant 
dependence on the number of blocks in this region. 

There are two possible statements of the dimensional reduction about the 
Polyakov loops 


Weak: the lowest state coupled to L(x ) is the same in both dimensions, 
so that the two functions must have the same shape in r at large r 
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Figure 2.1: Polyakov loop correlations Ld(t), D — 2 + 1 (circles) and 2 
(squares), for /? 3 equal to 29 (left) and 84 (right). Distance r is given in lattice 
units. In both cases L s = 32 and L 0 = 4. The 2D data were produced using h 2 
from Eq. ( |2.14| ). The continuous and the dashed lines result respectively from 
fitting formulae ( j2.34|) and (|2.33|) to the 2D data. 


• Strong: to the extent that the weight associated with is a good 
approximation to the integral over the non-static held of the (2 + 1 )D 
held (in the small p/T regime at least) the averages with the respect to 
the two weights of a static operator, such as L(r), are equal at large r 

We will show that the latter situation is approached at sufficiently high 
temperature. Using L 0 = 4 hxed, we investigate the temperature dependence 
of the correlations by varying /? 3 in the range [21,173]. High temperature 
means [3 3 sufficiently larger than the transition point in (2 + 1 )D. These /3 3 
values correspond roughly to temperatures in the range T/T c G [1.43,11.73] 

On the Fig.|2T| the correlations L 2+1 (r) and L 2 (r ) for /3 3 = 29.0 and 84.0 
are presented. We see that the two correlations are extremely close to each 
other: not only they have the same shape at large r, which is the primary 
prediction of dimensional reduction, but also nearly the same normalization. 
As announced, this finding favours dimensional reduction in the strong sense. 

In addition, we notice that this agreement between the two ways of com¬ 
puting the Polyakov loop correlation extends down to fairly small values of r. 
Since the normalization is set by | (L(x)) | 2 , it shows that, although it is lo- 
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cal, the Polyakov operator is not very sensitive to the short wave length terms 
omitted in the effective action. In Table |2.1| we present the results for (L 2 ) 
and (L 2+ i) to demonstrate this property. Relatively to their distance to one 
(their common value at f3 3 = oo), their difference, already ~ 6% at /3 3 =21 
(T/T c = 1.43), decreases to only about 0.6% at (3 3 = 173 (T/T c = 11.73). 


ft 

1 (l 2 ) | 

(£ 2 + 1 ) 

21.0 

0.56002 (62) 

0.53467 (16) 

29.0 

0.67007 (25) 

0.66120 (13) 

42.0 

0.76397 (20) 

0.76130 (12) 

84.0 

0.87392 (11) 

0.87435 (6) 

173.0 

0.93494 (13) 

0.93530 (11) 


Table 2.1: The average Polyakov loop as a function of /3 3 in (2 + 1 ) D and 
2D. 

In Fig. [2.2| , we illustrate more in details, for the (5 3 values not reported on 
in Fig. p.l|, how well the 2D and (2 + 1 )D correlations compare. Their ratio 
L 2 +i/L 2 is plotted against the distance r in lattice units. These data definitely 
support the statement that L 2+ i/L 2 remains quite flat and close to one for all 
distances and temperatures. Recalling that the lowest T value is only ~ 1.5 
times T C) and that the correlation functions decrease in r by about three orders 
of magnitude, we conclude that the effective local 2D action reproduces the 
(2 + 1 )D Polyakov loop correlations with a remarkable accuracy, soon above 
the transition and down to distances even shorter than 1/T. 

2.4.2 Scaling test 

As for every lattice calculation, the question arises if the lattice size is sufficient 
for the study, or finite size effects influence the measurements. This is normally 
checked by running the same simulation on the larger lattices and comparing 
results. The parameters of the bigger lattice are defined as follows. Given the 
temperature, the continuum limit is approached by taking 

r fixed L 0 = large (2.32) 

L 0 g* aT 

Scaling region is then defined as the region where physics of the model at 
fixed r does not depend on L 0 , if To is large enough. To check if we are in 
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Figure 2.3: Comparison of the Polyakov loop correlations in the 2D model 
for the two sets of parameters corresponding to T/T c = 1.97 [/3 3 ,L 0 ,L S ] = 

[29,4,32] and [58,8,64] i.e. for constant r Q2.20| ). The values of r for L 0 = 8 
(circles) are scaled down by a factor two in order to maintain the same physical 
scale. 


23 













Figure 2.4: The ratio between the two Polyakov loop correlation functions 
presented in figure [2~4| . 


this region, we compared the Polyakov loop correlations in 2 D for two sets of 
the lattice parameters, namely [/? 3 , L 0 , L s \ = [29,4, 32] and [58, 8, 64], Given T, 
doubling L 0 means dividing the lattice spacing a by two, so that the physical 
spatial size aL s of the lattice is unchanged. The resulting correlations are 
presented in Fig. |2.3| showing a very similar shape as a function of the physical 
distance. It is again instructive to consider their ratio, found to be quite flat at 
all distances (Fig. [2.4| ): within errors, there is no sizable deviation from scaling. 

The other observation is that this constant ratio is not one. And it should 
not be because it is given by the corresponding values of | (L(x)) | 2 . They are not 
physical quantities due to explicit (logarithmic) dependence of the 2D effective 
Lagrangian on a via its counterterm Eq.([2.28|). 


Hence scaling is verified in the range of interest. It justifies keeping Lq = 4, 
which is less expensive in computer time and allowed us to use the existing 
data of |28| in (2 + 1 )D. 
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2.4.3 Screening Lengths 


The fast decay of the correlation functions with r is consistent with the ex¬ 
istence of a finite spatial correlation length, fs, associated with the quantum 
numbers of the Polyakov loop operator, a color singlet scalar. A mass can be 
defined as usual by M s = ^g 1 . 

For the (3+1) dimensional model D’Hoker advocated the scenario (based on 
perturbation theory) that the screening mass Ms is 2M E , twice the so-called 
electric screening mass, because the lowest “state” coupled to the loop is a 
two electric gluon state j^Tj. In this case the L(r) correlation is expected to 


be proportional at large r to the square of the correlation function for a one 
particle state of mass M E : 


- (2m E ) 
J D 


(r) ~ d T 


[m E r] 1 / 2 


D -m E r 


+ 


[m E (L s - r)] 1 / 2 


0 -m E (L s -r) 


(2.33) 


In the present case, however, the infrared sector is more complicated and 
perturbation theory is doubtful. One of the options is that there exists in 
the Polyakov loop channel an independent screening mass Ms, associated with 
a true excitation of the 2 D system. In this case we parametrize the data 
according to 


L °‘” W “ c (p+F e "” Sr + ' (2 ' 34) 

The second term in both cases arises from the finite size of the lattice and its 
form is determined by the periodicity. It actually does improve fitting process 
significantly. The m symbols denote masses in lattice units, i.e. m = aM. 
These parametrizations respect the lattice symmetry r —> L s — r. Since we 
have no direct access to m E , the two expressions above differ in shape through 
the prefactors, r -1 / 2 and r~ l respectively. We checked that for r > 1 lattice 
artefacts are negligible in the mass range considered. This was achieved by 
comparing Eq. ( |2.34j ) to the lattice propagator 

LLatt(m s ,r) = cos (Pi r )LLatt(rn s ,p), (2.35) 

S P1,P2 

L~ L l att {rns,P) = f + 4sinh(m|/4). (2.36) 


We carefully analyzed our numerical data, and we find that the ansatz (|2.34| ) 
is by far the best one, giving a good fit of the data down to small r values. For 
an illustration, the continuous curves of Fig. |2.1| are fits of the form ( 2.34 ) to 
L 2 (r),r > r min , with r min = 4, and changing r min from 3 to 6 does not change 
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the resulting ms within errors. At the shorter distances a noticeable deviation 
exists, which can be either due to the presence of the more massive contribution 
or just be an artefact of the lattice UV regularization. It is anyway too small 
to be reliably analyzed. On the contrary fits to the same data of the form 
Eq. fl2.33D lead in Fig.|d]to the dashed curves, which are clearly not acceptable. 
The results of systematic fits of expression (|2.34| ) to all available correlations 
L 2 (r) and Z/( 2+1 )(r) are presented in Fig|h5|. We define the corresponding 
dimensionless quantities 


M s /(g 3 VT) = m s \jL 0 /3 3 /6 
gl/T = 6L 0 /p 3 


(2.37) 

(2.38) 


and plot them one versus the other. 


It is also interesting to check to which extent the calculated coefficients are 
important for the simulations and if they have only minor effect. Thus we 
also plot the masses resulting from a few simulations made with the classically 
reduced effective action (h 2 — — 0). They are far away from the former, 

showing the important effect of taking into account the non-static degrees of 
freedom. One may also want to check if the simulation which uses the so-called 
scaling form (|2.14|) of h 2 differs much from the simulation with the value ob¬ 
tained numerically from Eqs. (|2.18| , |A.52| ). The results are similar, especially 
at high temperature, but anyway distinguishable within our statistical accu¬ 
racy. This shows a great sensitivity of the static properties to the quadratic 
counterterm (|2.28|), an interesting feature per se. 


The scaling properties observed in the previous subsection reflect them¬ 
selves of course in the mass values. For the two cases compared there, we 
hnd M 5 a=0.331(6) and 0.169(3): the lattice spacing is reduced exactly by 
two within errors. Of course the fact that dimensional reduction works well 
for masses directly follows from its success for correlations. We preferred to 
illustrate it Erst on directly measured quantities, as done in Fig.(|2.3|). This 
comparison does not depend on any interpretation of the nature of the ob¬ 
served screening lengths. 
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Figure 2.5: Physical screening masses Ms in units of g%\/T versus g\/T , in (2 + 
1 )D (black points) and 2D (squares). Also shown are the masses obtained with 
the numerical value of h -2 from Eq. ( |2.18|) (triangles) instead of its asymptotic 
expression (|2.14|) , and with the classically reduced (labeled naive) action: no 
Higgs potential, /i 2 = ^4 = 0. Masses of the original model vanish at critical 
point, while in reduced one they do not 
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2.5 Spatial Wilson loop 

The spatial Wilson loops are measured on a lattice via 


W(R, R', x, y, m) = JJ U(x + i, y, 0) m JJ U(x + R, y + j, l) r 


(2.39) 


II U^x + i,y + R', 0) m J] rf(x,y + j,i) m (2.40) 




^ L s — 1 




(2.41) 


m= 1 £,y=0 


where N m is the number of measurements, and index m refers to the particular 
measurement. 

In our simulations the Wilson loops were measured every ten sweeps and 
blocked by 50. Integral autocorrelation time for the measurements of 14 x 14 
Wilson loops was, in the worst case, of the order of one and was negligible in the 
majority of cases. The measurement routine is quite expensive and consumed 
almost half of the CPU time. 

In order to extract the string tension, local potentials were first extracted 
from the ratios : 

T/Y (W(R,R')) 


V(R, R') = log 


(W(R, R' + l)) 


By definition, the potential is 


V(R)= lim V(R,R'). 

R '—>oo 


(2.42) 


(2.43) 


For each given R, V (R, R') was found to decay exponentially to a constant in 
R’. In practice, the constant regime is reached within errors above R' = 3 and 
the potential was fitted in the range R' G [4,12]. The string tension o was then 
obtained from the ansatz 

V (R) = V 0 + aR. (2.44) 

The errors on the string tension were calculated using the triple jack-knife 
method as follows: 


1. split the data, into N blocks 

2. omitting one block at time, produce N sets, each consisting of (N — 1) 
blocks 
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Figure 2.6: Potential extracted from the spatial Wilson loops, = 42 


3. in each of the sets repeat this procedure, so that we get TV x (TV — 1) sets 
each consisting of (TV — 2) blocks 

4. use these in traditional jack-knife way to obtain TV x (TV — 1) sets of 
V ( R , R') with errors 

5. repeat it to obtain TV copies of the V (. R ) with errors (which we need to 
do the weighted fitting) 

6. fit to the formula (|2.44|) get string tension and errors 


A typical form of the potential is presented in Fig. |2.6| . A fit using the 
linear potential (|2.44|) is sufficient to obtain a stable value of a within errors. 
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Figure 2.7: We plotted string tension for 4 different upper cutoff values against 
lower cutoff value 


To prove that, we produced for every /3 3 a fit from R m i n to R max ■ The R m ax 
values, depending upon the value of /3 3 , varied from 10 to 14, the R m i n varied 
from 1 to 4. To judge about stability of the fit we plotted the resulting value 
of cr with errors as the function of R , min , R ma x■ As it may be seen from Fig.[2vj 
the fit is indeed stable. This is in good agreement with fluctuating flux tube 

0 


model 


the Luescher 1/R term is not expected to be present above the 
phase transition. 

The x 2 value in the 3D case was of the order of one and about ten times 
less in the 2D case, likely due to the large correlation between points. 

For the dimensionally reduced treatment the case of the spatial Wilson 
loop W in 3D is more difficult than the Polyakov loop correlations. It is not a 
static operator, and thus by no means is predicted to be reproduced well in the 
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Figure 2.8: The square root of the physical string tension in units of g^VT as a 
function of g\/T in (2+l)D (filled circles) and 2D (squares). The line denotes 
the scaling limit a/2/3 in 2D pure gauge theory (|1.21|) . 


dimensionally reduced model in the high temperature region. However, there 
are the following considerations: 


i) a non-zero spatial string tension <j- 2 +i is known to exist |28| above the 
deconfinement temperature 


ii) the pure gauge 2D theory is confining and produces a finite string tension 
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iii) in two dimensions we expect confinement and thus a finite string tension 
<j 2 to survive when the Higgs field (j) is turned on. 


It is thus interesting to compare these three quantities as a function of the 
temperature of the (2 + 1 )D model. For this comparison we take a?cr 2+1 and 


31 
























a 2 cr 2 from the analysis of W data provided by |2B[ and by our simulation, 
whereas a 2 ^ is computed analytically in Appendix 2: 


a&1 “ ft + M + 0(/?2 
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The next terms in this expansion are easy to derive using the computer algebra 
system. However, they are not required for our purpose. This expression can 
be equivalently rewritten 


an 2 7 g\ . 2 

2 = - + — %{aT) 2 + O 


glT 
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showing that up to scaling violations of order 

(' aT? = T (2.45) 

L, o 

<7° scales as 2g|T/3 at finite T . The quantities reported in Fig. [27^ versus gl/T 
are the values of y/a in units of g^VT, which we thus compare to 2/3. The 
numerical values of cr|/(.g^T) and og+i / [g 2 T) were obtained from the simula¬ 
tions for Lq = 4. 

All three string tensions are quite close to each other; they differ by at 
most 10% within our temperature range. From the two-dimensional point of 
view, we see that introducing the 0 held modifies the picture quite weakly: 
no sizable g 2 /T dependence around a value close to ^/2/3. The behavior for 
a 2 + i is different, with a sizable slope in g\fT. The difference observed between 
a 2+ i and a% may be understood as a consequence of W not being a static 
operator: the average of its non-static modes with the (2 + 1 )D weight are 
missing in its calculation with the effective action. However, scaling violations 
may also contribute differently to these two cr’s, so that def ini te conclusions 
require complementary simulations. 

Another possible reason for this disagreement are the heavy space correla¬ 
tions in the Wilson loops measurements, which may lead to systematic errors 
for the a measurements. Our data, presented in Fig.[2~E|. do not lead to the 
definite conclusion whether or not the infinite T limit of the string tension in 
the two- and three-dimensional models is the same, and if it is equal exactly 
to 2/3. 
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Chapter 3 

Reduced Model Per Se 


3.1 Phase structure of the model 

As we have seen in the previous chapter, the reduced model reproduces the 
(2 + 1 )D model very well; it is interesting thus to study the properties of 
the reduced model as such. As we pointed out in Chapter 2 , the wanted 
deconfinement phase transition is absent in this model; one can see it either 
from the fact that screening masses do not vanish, or via direct test. The reason 
for that is the explicitly broken Z( 3) symmetry in the perturbative approach. 
One of the remaining symmetries is the R r (see |2.30|) - the imaginary time 
reversal symmetry of the original model. As it was pointed out in the studies 
of the (3 + l)H case, there exist two phases, symmetric and broken with respect 
to this symmetry. The same scenario happens also in our model. The test runs 
on the small lattices have shown the possibility of the very different dynamics 
of the system. Indeed, we were able to see a strong first order phase transition 
happening. The order parameter for this transition is 

0\ = {trA 3 ) (3.1) 

It is equal to zero in all the simulations discussed previously and is different 
from zero in the broken phase. 

We performed extensive simulations for lattice sizes 8 x 8 and 6 x 6 with 
parameter L 0 = 4. To see how the transition point moves with scaling we also 
simulated on 8 x 8 , 10 x 10 and 12 x 12 lattices with L 0 = 8 . Couplings /3 2 and 
hi were fixed at three-dimensional values corresponding to (3 3 = 29 (/? 2 = 116, 
hi = 0.954930) and = 58 (/? 2 = 464, hi = 3.819718) respectively. The oper¬ 
ators of interest this time were TrA 2 and TrA 3 . To map the transition point 
we performed a large number of measurements, which varied from 3.000.000 to 
30.000.000 depending upon coupling. The typical autocorrelation time was of 
the order of several thousands due to closeness to the phase transition. 
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We used the histogram technique to find out the position of the phase 
transition. Unlike in the previous chapter, the data were not blocked to increase 
the statistics for histogramming. Typical histograms for TrA 2 are presented 
on Fig. |3.1| and Fig. |3.2| for L 0 =4 and 8 respectively. They have a clear two- 
peak structure, which is the traditional picture for the phase transition of the 
first order. Behavior of the system was significantly different in these cases: for 
L 0 = 4 already at 6 x 6 lattice a very sharp two-peak structure was present, and 
two phases are separated clearly. For Lq = 8 the transition is much softer, the 
relative height of the peaks is smaller compared to the clip. This is compatible 
with the fact that we decreased the lattice spacing twice when increasing L 0 . 



Figure 3.1: Histogram, TrA 2 ,6 x 6 lattice, Lq = 4 

A typical histogram for ()\ is presented on Fig. |8.-'i| - a clear three-peak 
structure is present; the scaling behavior also turned out to be similar. This 
histogram shows that in fact there are not two but three different phases, 
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0.0025 



Figure 3.2: Histogram, TrA 2 , 8x8 lattice, L 0 = 8 


the symmetric one and two broken phases with positive and negative value of 
the order parameter. For the data analysis we used the Ferrenberg-Swendsen 
multiple histogram reweighting algorithm. It allows to use the data from several 
histograms, effectively accumulating statistics for a more precise determination 
of the critical point. Details of the algorithm may be found in [[IT]]. The single 
histogram reweighting method was used for the cross-checking. 


There are several criteria which one may use in case of the first order phase 
transition, among them the equal weights criteria and the susceptibility criteria. 
To employ the equal weights criteria we need to take into account that there 
are three phases. Thus the criteria should be formulated as such that broken 
phases weight is 2/3. After scanning the h 2 space, we spotted the approximate 
point where transition happens and performed about a dozen of long runs 
in that area. Then we used the Ferrenberg-Swendsen reweighting to produce 
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Figure 3.3: Histogram, TrA 3 , 8x8 lattice, L 0 = 4 


the table of the pairs ( h 2 ,Wb ) where H/, is the weight of the broken phase 
(i.e. number of measured values which are on the left side of the minima, 
divided by the overall number of measurements) with a very small step in h 2 . 
Afterwards the hi corresponding to the closest to 2/3 value of Wb was selected 
to be an estimator for the critical point. This estimator was used in jack-knife 
like treatment to obtain the errors on the critical value of hi. 

A typical picture for the susceptibility is presented on Fig.[T4] where sus¬ 
ceptibility is plotted versus h 2 for L 0 = 4. To estimate the h critica i, we find the 
absolute maximum for every jack-knife bin, ignoring the actual errors on sus¬ 
ceptibility. As it can be seen from table |3.1| below, the susceptibility maximum 
criteria estimate gets closer to the equal weights criteria one when increas¬ 
ing volume and, moreover, coincides within errors for the bigger value of the 
parameter L 0 . 
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6x6 

8x8 

oo 

3D 


equal weights 

-5.186(2) 

-5.225(4) 

-5.275(12) 

-6.071 

susceptibility 

-5.165(3) 

-5.208(4) 

-5.263(13) 

-6.071 


8x8 

10 x 10 

12 x 12 

oo 

3D 

equal weights 

-7.28(2) 

-7.33(2) 

-7.37(2) 

-7.45(9) 

-8.06 

susceptibility 

-7.28(2) 

-7.34(2) 

-7.39(4) 

-7.44(11) 

-8.06 


Table 3.1: Critical couplings for L 0 = 4 and L 0 = 8 


We also provide in the table the reduced values (labeled 3D) and two ex¬ 
trapolated values, an “infinite volume” value and the 32 x 32 value, which we 
need for making the conclusion about where we are on the actual lattice we 
work on. The scaling behavior of the following type was assumed 

h v = h in{ + + O (3-2) 


To check the validity of this ansatz we compared the predicted from the two 
points value for the third point with the actual value, and they agree within 
errors. 

The phase diagram is sketched on Fig.|3.5|, the dotted line is there to guide 
the eye. The reduced values are obviously in the broken phase, and situation 
does not get any better with increasing L 0 . The Polyakov loops and the string 
tension are not reproduced by the simulations in that phase, ffowever, due 
to it being a very strong first order phase transition - it never (on a scale of 
several gigaflops-monthes) happens on the 32 x 32 lattice which we used in 
the previous chapter. Starting from a weak held configuration, we are in the 
symmetric phase and always stay in it. 

I3| as done 


One could perform here a more detailed analysis following Ref. 


in the determination of the critical exponents and transition temperature in 
|I2| along with the study of the larger lattices. Latter would require some 
advanced simulation method (e.g. multi-canonical simulations), ffowever, as 
long as the transition is of a strong first order we conclude that for the small 
values of L 0 our lattices are big enough to make a statement that the reduced 
values are in the broken phase. 


3.2 The spectrum of the model 

In this section we will study the lowest states of the spectrum of the reduced 
model per se. We restrict ourselves to the operators directly related to the 


38 





























Figure 3.5: Phase Diagram of the Model 
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static operator L{x) because only they are directly relevant for the (2 + 1 )D 
system. 

Let us now expand the L{x) in powers of A{x) and study the operators 
A n (x ) and connected correlations of their traces. We dehne 


A n (x) = Tr A n (x), 

A nt m(x) = {A n (x) A m {Qt)) - (A n )(A m ). 
A n = {A n {;x)) 


(3.3) 

(3.4) 

(3.5) 


Any operator A n {x) is gauge invariant, even or odd under the /?, r -symmetry of 
the 2D action for n even or odd respectively. We will continue working in the 
unbroken R r phase only, where A 2p + 1 = 0. 

Using the same action as in Chapter 2 we have performed a numerical 
simulation of the two dimensional model for lattice size L$ = 32 till the end 
of the chapter. The (3 3 values are 29, 42, 84 and 173, while L 0 = 4. They 
correspond to the values of T/T c equal approximately to 1.97, 2.85, 5.70 and 
11.73 respectively. We were able to extract information from the operators and 
correlations corresponding to n = 2 to 5. The cases n — 2 and 3 for the 3 D 
reduced model were investigated in p0 |. 

Here we will present our results for the A n ^ m correlations measured and 
describe them for each temperature in terms of two states, which we will call 
S and P, respectively even and odd under R T . They appear for n, m both even 
and both odd. Their physical masses will be denoted Ms and Mp. 

We will use in this and the next section either the bare lattice parameters 
r and /? 3 or the physical quantities R and T, which are related to each other 
by the following equations: 


RT = r/L$ (3.6) 

We show in Fig.[L6] the on-axis correlations A n>m (r), n < m G [2, 5] at 
T/T c = 1.97 (/?3 = 29). They are plotted against RT, that is the physical 
distance in units of the inverse temperature. As expected from general con¬ 
siderations one can see two types of behavior, for the even and odd channels. 
I 11 the even cases, the three correlations all have the same shape, and they 
decrease by about one order of magnitude each time two more powers of A are 
involved. The same is true for the odd cases, with a common decay of the cor¬ 
relations steeper than in the former case (smaller correlation length in lattice 
units). The overall situation is similar for /? 3 higher. For n,m larger than 5, 
as well as for /? 3 very large, the signal/noise ratio becomes very small. This 
can be understood at the qualitative level by noting that the rescaling Eq. (p.8|) 
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Figure 3.6: The on-axis correlations A n ^ m {r ) at T/T c = 1.97 (/5 3 = 29), versus 
the distance in units of 1/T. The even cases [n,m] = [2,2], [2,4] and [4,4] all 
have the same shape, and the odd cases [3,3], [3,5] and [5,5], again similar 
with each other in shape, are steeper. 
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of the H-helds normalizes the kinetic term for the 0-helds to the standard, 
parameter independent form 1/2 Tr{D i <^)‘ 2 . Hence, if the held renormalization 
by the interactions is weak, the (j) correlations should depend only weakly on 
/3 3 , which means that A n ^ m scales like f5^ n+m ^ 2 . This will be illustrated more 
quantitatively in the next subsection. Due to this scale factor, the H-helds 
remain “small” in practice down to quite low values of d 3 , which a posteriori 
explains why the perturbative reduction may still work at a temperature as low 
as 1.5 T c . In fact, we checked that the Polyakov loop correlations are actually 
fully reconstructed within errors by keeping {n,m} up to {5,5} only in their 
expansion in H n m ’s obtained from the small A expansion of ( |2.31|) . 

Now we want to analyze these m data quantitatively in terms of the 
lowest states of the spectrum. Let m, = a M t be the lowest mass with quantum 
number i = S, P, in lattice units. For particle i with momentum p we introduce 
a lattice propagator A Latt (rrii, p) in momentum space: 


A Lait("h>.P) = f + 4sinh(m 2 /4), 

p 2 = 4sin 2 (pi/2) + 4sin 2 (p 2 /2). 

The corresponding contribution to A njTn (r) then reads 

An,m{rrii,r) = g l n , m -j^ ^ cos(pir)A Latt (mi,p), 

S Pl,P2 


(3.7) 


(3.8) 


where g l nrn measures the residue of A ntfn at the pole of (|3.7|), which on large 
enough lattices sits at p 2 ~ p 2 ~ —mf. With our definitions, g l nm is non-zero 
only for i — S if n and m are even, and for i = P if n and m are odd. 

There is a simple result which one may get directly from Eq. (|3.8|) . For the 
case of the single-particle propagation between two spinless states the residue 
g l nm factorizes (see, for example, [[57]]): 




= 7n7r, 


(3.9) 


This property we can probe directly on the correlations since, as r becomes 
large, it implies 


_ A nt n{T^ H n _|_2,n+2 (j’’) _ 

Al tn+2 ( r ) 


(3.10) 


We demonstrate this fact for the temperatures T/T c = 1.97 and 5.7 (/? 3 = 29 
and 84) on Fig. [j.7| - {rTT]| .The displayed quantities A" 2 and X 3 (symbols o) indeed 
approach one for the large distances in all cases. However, the quality of the 
data is poorer for A" 3 due to the correlations involving A 5 getting very small. 
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Similar results are obtained for other values of T/T c . We thus conclude at 
this point that a single particle propagation accounts very well for the largest 
correlation length occurring in each of the two channels. Most of the observed 
deviations of X n from one will be interpreted in the next section in terms of 
the two particle state contributions (symbols o in the same figures). 

We now proceed to assign values to the two lowest masses Ms and Mp 
expected from the above findings. This we do by various ways in order to 
further enforce the statement that the correlations do have the characteristics 
associated with the pole structure of Eq.(|3.7|). Down to r ~ 1, an excellent 
approximation to the on-axis correlation ( |3.8|) is given by 


A, 


n,mym h r) ~ c 


[mjr ] 1 / 2 


e~ mi r + 


[m,i(L s - r)] 1 / 2 


0 —m.i(L a —r) 


(3.11) 


where c is constant in r. The second term takes into account the finite size of 
the lattice. 

We performed fits of this formula to all our A n ,m( r ) data taken at r > r min . 
These fits are stable with respect to r min , provided it is larger than about 4, 
and the values found for rri, in different correlations are always consistent with 
each other. The smallest errors were obtained by using fits to A ,2 and A 3 ,3 for 
the reasons described above. 

An alternative way for the extracting the effective masses without any fit¬ 
ting is the use of the 0—momentum correlations. They defined for a generic 
x —space correlation C(x 1 ,^ 2 ) by 


C°(r) = 


C( Ai x i)- 


X2 


If the lowest mass in C is m, the ansatz (|3.7| ) gives 

C°(r) oc cosh (m(L s / 2 — r)), 


and m can be extracted at any r by inverting this relation: 


m — log (Vl'r i + yT-ir] - 1 j , 

, C°(r + 1) + C°(r - 1) 

r) = - 2 CW) -' 


(3.12) 


(3.13) 


As an overall consistency check, we have extracted an effective mass m e ^ (r) 
from 0—momentum Polyakov loops correlations ( P.31Q and compared it to the 
ms values obtained by our fits to A 2 , 2 - We find that 7n e ^(r) is indeed nearly 
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T/T c 

Ms/ValT 

Mp/^gjT 

1.97 

1.39(3) 

2.51(13) 

2.85 

1.48(5) 

3.02(16) 

5.70 

1.89(5) 

3.22(15) 

11.73 

2.12(8) 

3.44(11) 


Table 3.2: Masses in units of a JgfT for the S and P states, as measured from 
fits to A 22 and A 33 respectively, for different values of T/T c . 


constant, in fact slowly decreasing towards a value compatible with ms , due 
to smaller and steeper contributions to ( j2.31| ) of the heavier particle P. 

We also investigated the D’Hoker’s perturbative scenario, as we have done 
in Section |2.4.3| . Again the largest correlation length in A n n should be n times 
shorter than the “Debye screening length”, the inverse of a mass m E associated 
with “electric” gluons of the initial (2+l)D model (the scalars of the reduced 
model). If such was the case, the on-axis correlations A n>n should rather look 
like 


A n ,n(nm E ,r) oc 


[m E r] 1 / 2 


D -m E r 


+ 


[m E (L s - r)] 1 / 2 


-,-m E (L a -r) 


(3.14) 


which differs in shape from (|3 .1 1|) , as was illustrated in the previous chapter 
for the Polyakov loop correlations. We, nevertheless, tried fits with ( |3.14 ), but 
got a definitely worse agreement in the range of temperatures, which we have 
investigated, i.e. up to 12T C . Hence, this scenario is ruled out by the data in 
this temperature range; if a mass can be defined for the electric gluon in high 
temperature QCD 3 , it is most probably larger than both ms /2 and mpj 3. In 
the ’’constituent gluon” picture, as advocated in Ref. f58|, one would have 
bound states instead of a cut. One would, however, expect mp/ms ~ 3/2. 


Our final results for the S and P masses in units of the scale \fgfT are 
collected in Table [3]2] for the values of T/T c investigated. They are taken from 
fits to A 2) 2 and A 33 respectively. The values for M s agree with those obtained 
before from the Polyakov loop correlations. As can be seen from the tables, 
the ratios Mp/Ms vary with T/T c . There is, however, no clear tendency in 
the region we have investigated, the ratios being 1.8, 2.0, 1.7, 1.6 in order of 
increasing temperature. This ratio may, of course, go to 1.5 at still higher 
temperatures. 
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3.3 Weak “Strong” Interactions Between Col¬ 
orless States 

Here we will show that even at quite short distances (r small compared to rrig 1 ) 
all the condensates A n = (A n (x)) and correlations A nim (r) can be reconstructed 
to a good accuracy from the data for A 2 , A 2i 2 and A 3i3 . The assumption is 
that the elementary fields A a (x ) (Greek superscripts are color indices) interact 
only through S and P exchanges between the non-interacting composite A 2 (x) 
and A 3 (x), the scale of the fields being fixed by the size of A 2 , while A 3 = 0. 
The details of the derivations are presented in Appendix 3, while here we limit 
ourselves to the simplest applications and give the results, starting with the 
local condensates. 


3.3.1 Weak Residual Interactions: The A-fields conden¬ 
sates 

Since SU(3) has rank 2, any A n (x) can be reduced to a polynomial in A 2 (x) 
and A 3 (x). For n odd A n is zero by R r symmetry. For n even, we apply Wick 
contraction to all pairs of A a elementary fields, followed by the meanfield-like 
substitution 

A a (x)A f 3 (x) -> \s at pA 2 (x). (3.15) 


As an illustration consider A 4 . With the definitions of Section 1373, we have 


2 A 4 (x) = Al(x) = — ^ A a (x)A a (x)A f 3 (x)A , 3 (x). (3.16) 


q,/3=1 


There we apply ( [3. 15 ) and then replace A 2 (x) by its average A 2 . The A a A a 
and A 13 A 13 contractions give (8 x A 2 / 4 ) 2 , and the additional contributions from 
a = P give 2 x 8(A 2 /4) 2 . Noting that 


:(0) = (Al(x)) - A\ 

(3.17) 

be put into two equivalent forms 


2A 4 = - A 2 , 

(3.18) 

^ 2 , 2 ( 0 ) = - A\. 

(3.19) 


This prediction is remarkably well verified in all cases. At p 3 = 29, the left 
and right hand sides of ( f>.19| ) are respectively 4.86(1) 1CT 3 and 4.825(10) 10 -3 . 
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They are 6.093(13) 10 4 and 6.069(1) 10 4 at fa = 84. The formula for the 
A 3)3 (0) correlations 

-43,3(0) = (3,20) 

leads to the similar correspondence. In this case the left and right hand sides 
are measured to be 2.105(7) 10” 4 and 2.095(6) 1(T 4 for fa = 29, 9.365(30) 10~ 6 
and 9.344(3) 10~ 6 for fa = 84. Hence, the effects of residual interactions via 
higher orders effective couplings in A 2 (x) and A$(x) are less than the percent 
in the correlations at zero distance. 

Before going to the correlations at non-zero distance, let us discuss their 
normalization, as measured by the values of A 2: 2 ( 0 ) and A 3 i 3 (0) described just 
above. At the beginning of the previous section, we argued that the behavior 
in fa, n , m observed for A n)m could follow from the absence of a large renor¬ 
malization of the 0-fields (defined by Eq. (|2.8|) ) by the interactions . Here we 
note that in the confined phase the effective degrees of freedom are the massive 
composites 0, = Tr0*, i — 2, 3, so that in the limit where they are considered 
as free fields, one may write (see ® 3 » 


(A(o)A(o)) - Ri / <*V 


p 1 + 4 sinh(mf /4) ’ 
p 2 = 4 sin 2 (pi/2) + 4sin 2 (p 2 /2), 


(3.21) 


the residue Ri being one if neither 0 nor the composites get renormalized. 
We computed Ri as the ratio of the l.li.s. of (|3.21|) , directly measured, to 
the integral in the r.li.s, evaluated numerically on the lattice for the mass 
values fitted to the correlation data. The result is shown in Fig. |3. 1 1| : in the 
whole temperature range, both residues in the even and odd channels remain 
uniformly very close to one. 


3.3.2 Weak Residual Interactions: 
Correlations 


Properties of the 


As we have seen in section |3]^ (see Fig. |3~2| ), the different A nir7l (r)’s correspond¬ 
ing to the same channel have very similar shapes. Their analysis in terms of 
one particle exchange was successful, confirmed by the agreement with residue 


factorization. Nevertheless, although the quantity X n of Eq. (|3.10|) does go to 
one at large distances, it is significantly different from one at medium and short 
distances (see Figs. |T7|j8. 10| ). We will now show that two particle exchange is 
responsible for most of this lack of factorization. 

The simplest consequence of our assumptions for correlations at finite r is, 
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Figure 3.7: Residue factorization: data for the quantity X 2 , Eq. ( |3.10| ) at 
T/T c = 1.97 (/? 3 = 29) versus the distance in units of 1/T. It approaches one 
at large distances. The quantity X 2 corresponds to our interpretation (Eq. 
( |3.24| )) of the deviation from one of X 2 at shorter distances. 
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Figure 3.8: Residue factorization: data for the quantity A 3 , Eq. ( |3.10| ) T/T c = 
1.97 (/3 3 = 29) versus the distance in units of 1/T. ft approaches one at large 
distances. The quantity A" 3 corresponds to our interpretation (Section 4, Eq. 
(|3.25|)) of the deviation from one of A" 3 at shorter distances. 
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Figure 3.9: Same as in Fig. ^77] at T/T c = 5.7 (/3 3 = 84). 
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Figure 3.10: Same as in Fig. |3]^ at TjT c = 5.7 (/3 3 = 84). 
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using Eq. ( |3 .1 8| ) , 



(3.22) 


which is very well verified at any distance as shown in Fig. |3.12| for T/T c = 1.97 
(/?3 = 29.0). A similar agreement is found for the relation 


^ 3 , 5 ( r ) — 777 A 2 A 3 i 3 (r), 


(3.23) 


derived in the appendix. A new situation arises when we consider A 4j4> or A 5i5 
where both the initial and final states may couple to a two particle state, (S S ) 
or (S P) respectively. Then the intermediate state in a connected correlation 
between 0 and r may consist of either one or two particles. For example, to 
compute A 4 j 4 ! (r), we apply the substitution rule (|3.15|) to the sum ( |3.16|) , and 
then average using the definitions of A 2)2 and A 2 . One finds 



and 



These both expressions contain two parts. The first term in the square brackets 
is the contribution of the single particle propagation. The second term is a 
product of two propagators in space, and is that of a two-particle intermediate 
state. It provides a correction to exact factorization. 

From the definitions Eq. (|3.10|) of X 2 and A " 3 one actually gets the following 


estimates: 



(3.24) 


(3.25) 


The estimates X 2 (r), X 3 (r) are displayed in Figs. |3uj - |3Tl (resp. |T9| - |3.10|) , for 
comparison with the measured values X 2 (r), X 3 (r) at T/T c = 1.97 (resp. 5.7), 
i.e. f3 3 = 29.0 (resp. 84.0). We see that the corrections to factorization 
implied by the two-particle propagation provide a reasonable explanation of 
the behavior observed for the X’s at intermediate and short distances. This 
is especially true in the case of A 2 , showing that there is very little room for 
contributions from direct non-quadratic couplings in A 2 (tc) in the full effective 
action (that resulting from integration over the gauge fields). This justifies our 
statement that the residual interactions between the colorless boundstates of 
the adjoint scalars are very weak. 
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Figure 3.11: The residues Ri defined by Eq. 
whole temperature range. 


(|3.21|) stay close to one in the 
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Figure 3.12: Plot of 4/L 24 /5y4 2J 4 22 at T/T c 
one if Eq. (|3.22|) exactly holds. 


1.97 (/?3 = 29). This quantity is 
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Chapter 4 

Z(3) symmetric action 


In this chapter we describe an attempt to overcome the absence of the decon¬ 
finement phase transition in the model. As we mentioned before, this phenom¬ 
ena is not possible due to the fact that we broke the Z( 3) symmetry explicitly 
when doing the perturbation expansion. One may try to restore this sym¬ 
metry by hand and derive an effective action possessing it. We will use Ao 
variables throughout this chapter, them being a more natural choice for our 
consideration. The starting point is then: 


Sgff (U, Aq) — S^r (U) + Si n t (U, Ao) + Spot (Ao) 


where S%y ( U ) is the standard Wilson action in two dimensions as defined in 


Eq.(|TA3D, 



X l 


Dj A 0 = U (x; i) A{x + ai) U (x;i) 1 — A (x) 


S pot (A 0 ) = ^ hoTrAl (x) + h 4 ( TrA\ (x)) 2 


X 


with couplings 




li — In L 0 + - In 2 - 1 

We may now formulate the basic requirements for this action: 


i) should be Z 3 symmetric 
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ii) should have also the basic gauge symmetry 

iii) should reproduce the results from the ’’old” model at the large value of 

Ps/Lq, because we need it still for describing the (2 + l)-dimensional 
model 

One would also like that such a model predicts ^-symmetry restoration 
(reconfinement) as Pq/Lq (temperature) is lowered. The Polyakov loop should 
be the order parameter. Close to one (up to a Z 3 phase exp(±27ri/3)) at 
high temperature, it is zero (on the lattice it means - becomes small) in the 
symmetric phase. 

It is obvious that there are infinitely many ways to create the desired action. 
The most reasonable for us seems the following form: 

Sz 3 = S 2 w (U) + S uy + S v (4.1) 

The plaquette action S^y ( U ) remains untouched. The kinetic term Sjj,v we 


rewrite using the Wilson Lines 

V(x) = U^°(x) (4.2) 

where Uq is the element of the 577(3) group, which we place on every site of 
our lattice. It relates to the A 0 held as usually via Uq = exp(L4 0 ). The simplest 
choice for the kinetic term is 

Suv ~ — V'' ( 1- ?R.TrU(x; n)V(x + ajl)U(x\ /j,)~ 1 V~ 1 (x) J (4.3) 

V 3 / 

There are two invariants for the SU (3) group, so one can choose any two 
independent quantities as the basic variables for writing up the potential. We 
made our choice in favour of the following quantities: 

0(x) = \TrV(x) 0 *{x ) = \TrV*(x) (4.4) 

o o 

that is the Polyakov loop and its conjugate. This choice of variables leads to 
the following considerations 

• gauge invariance implies that Sy depends only on the set of all pairs 
{ 0 (x), 0 *(x)} 

• by Z 3 invariance, the lowest degree local monomial in 0 , 0 * is O 2 (x) = 
0 (x) 0 *(x) (degree two). 
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next terms should be 0 3 (x ) and 0 * 3 (x ) with degree three, there is no 
reason for them to be absent. They are independent degrees of freedom 


• reality implies that at order 3 in O the local action depends only on 

0 3 (x) = ^(0 3 + 0* 3 ). 

and terms proportional to it are expected on very general grounds in the 
full effective Polyakov loop potential. 

The simplest form of such potential is then 

Sv = X ] ( A 2 \ 0 { x )\ 2 + A 3 Re(0 3 (x))) (4.5) 

X&A.2 

The corresponding couplings A 2 and A 3 have to be fixed by again matching 
the small A 0 expansion of (|4. 1 |) to Sa 0 • Although this proposal still contains a 
substantial amount of arbitrariness, the following may be stated: 


i) All desired symmetries are implemented, 

ii) All independent degrees of freedom are involved, 

iii) In the region of phase space where the O’s are close to one, the action 

looks like the old one. The hope is that it is the relevant region at large 
(3 3 /L 0 values, 

iv) lu the region of phase space where the O’s are close to zero, the action is 

truncated to the two lowest degree monomials in O, O*, that is two and 
three, not two and four. The hope is that it is the relevant region at low 
(3 3 /L 0 values. 


The other option would be to ignore the requirements ii) and iv) and use 
instead of 5P(0 3 (a;)) the term with 0 4 (x). This way it is also possible to match 
all the terms of the potential (see below). 

One wants to be working in the unbroken phase, and assume that trA 3 = 0. 
The obvious choice for the kinetic term is 

Si = - 7 ^- (1 - R eTr {U (x, -/i) O (x + /x) U f (x,/u) O* (x + /i)} 

VjLq L ^ 


57 


Expanding the action ( |4.5| ) one easily gets: 

AC _ 9 20/i — 1 
2 ““ 16 vr Ll 
3 12p — 1 
As “ 8 vr Ll 

Here is the table for the coefficients for f3 3 = 29: 


L 0 


A^ 

4 

0.4631026079 

-0.182256888 

8 

0.1545592536 

-0.06107766263 


The model described above represents a good candidate for the improved 
effective model, which can allow us to calculate the three-dimensional quantities 
closer to the phase transition. Our test runs indicate that this model does have 
the second order phase transition in /3 2 with Polyakov loop being an order 
parameter. However, it still requires a detailed study, which goes beyond the 
scope of the current work. 
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Chapter 5 

Conclusions and outlook 


This work is the detailed investigation of the dimensional reduction for SU(3) 
gauge theory in (2+l)-dimensions at high temperature. We have produced 
an effective reduced model for the static variables following (T^. [15|. |l(| . This 
means that the integration over the non-static modes has been performed per- 
turbatively, and the effective couplings have been kept to one loop order for 
the two point and four point functions. Higher order contributions, neglecting 
which was a part of our approximation, are unimportant at large distances and 
high temperatures. 

We investigated the validity of this approximation for the correlation func¬ 
tion between Polyakov loops and corresponding screening masses, as well as the 
spatial string tension. For the (2+l)-dimensional theory and the corresponding 
2-dimensional reduced adjoint Higgs-gauge model we were able to obtain very 
precise numerical results. Therefore, we could make a detailed comparison, 
including a scaling analysis, which strongly supports the assumption that our 
lattices correspond to a large enough value of the temporal extent L 0 . 

For the correlation between Polyakov loops the dimensional reduction works 
very well down to temperatures of the order of 1.5T C , where T c is the critical 
temperature of the deconfinement phase transition of the original model. This 
is true even at distances down to or below 1/T . And it does even better 
than in the (3+l)-dimensional case. The correlation between Polyakov loops 
is well described by a simple pole in momentum space, even at high tempera¬ 
tures. This situation is not manifested in (3+l)-dimensional theory where at 
high temperatures the data favour the exchange of two Debye screened gluons 
|43| , [29| . [53], [44|. The disagreement may be related to the stronger infrared 
divergences in (2+1) dimensions, which invalidate the perturbatively resumed 
Debye screening picture. A further investigation of this behavior in (3+1) 
dimensions would certainly be very interesting. 

The spatial Wilson loops in (2+1) dimensions do not correspond to a static 
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operator and is not reproduced at all by itself. However, for the string tension, 
which is extracted from Wilson loops of extent larger than 1/T, one can hope 
that the non-static corrections are small. Our comparison favours this situa¬ 
tion. In fact, the Higgs sector seems to have little influence on the value of 
the string tension; there is a fairly good agreement also with the analytically 
solvable pure two-dimensional gauge theory where confinement is given by the 
two-dimensional Coulomb potential. The differences between the string ten¬ 
sion in the three cases considered are of the order of a few percent. To make 
a definite statement if these differences are real continuum effects, one must 
further study finite size and scaling corrections. 

Our work shows that it may be possible to explain the non-perturbative 
features of (2+l)-dimensional QCD in the deconhned phase with a relatively 
simple two-dimensional model. The model possesses, however, a non-physical 
phase transition associated with the symmetry of the Higgs held sign reversal 
R r , which corresponds to the time reversal of the original model. The reduces 
values of the Higgs potential are in the non-physical broken phase. However, 
we are able to simulate with these values in the metastable phase, and never 
drop into the broken one on sufficiently large lattices due to this transition 
being strongly first order. 

We studied the reduced model in the region of the unbroken symmetry 
to identify two boundstates S and P, respectively, even and odd under R T 
and thus coupled to monomials, respectively, of degree 2 n and 2 n + 1 in the 
higgs held. The S signal coincides with that previously obtained from the 
Polyakov loops correlations, where, however, the P-state contributions could 
not be disentangled. These results came from the measurement of three even- 
even and three odd-odd distinct correlations, as functions of the on-axis lattice 
distance r. We very carefully analyzed their shape in r, with the result that 
in all cases the signal found was that expected from the occurrence of genuine 
poles in momentum space. Our data invalidated the D’Hoker’s scenario where 
the decay of such correlations with the distance reflects the propagation in 3 D 
of p = 2 n or 2 n + 1 “electric gluons”. In other words, the perturbative result 
that the correlation length equals to 1/p times the “Debye screening length” 
is inadequate in the present case. 

By comparing the size of the three different correlations measured for each 
of the S and P sectors, we were able to show that residue factorization holds, 
as expected on general grounds when one particle propagates between different 
states. The agreement with factorization was expectedly found to be particu¬ 
larly good at large distances. We were also able to demonstrate that deviations 
at shorter distances are to a large extent compatible with propagation of two 
particles, namely two S or S and P, respectively, in the S or P channel. We 
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thus make a conclusion that the scalar sector of the reduced model at large 
distance can be described by the two weakly interacting colorless states. 

The perturbative approach used to derive the reduced model explicitly 
broke the Z( 3) symmetry of the (2+1) dimensional system, thus we do not 
have the phase transition which would correspond to the deconfinement phase 
transition of the original model. We presented here a candidate for the effective 
model, based on the original spatial gauge variables and Wilson Lines. This 
model is likely not only to be able to reproduce the full model at high tem¬ 
perature but also exhibits a deconfinement-like phase transition with Polyakov 
loop being an order parameter. Careful study of this model is needed to make 
a definite statement. 

Farther interesting topics for investigations are of course more real-world 
models. The similarly detailed analysis for pure gauge QCD in (3 + 1)D would 
be interesting, with extending to the full QCD with dynamic fermions. Fur¬ 
thermore, the study of thermodynamical quantities like pressure and entropy 
may be of interest, as well as study of the operators which are even and odd 
under other symmetries of the model - space reversal and charge conjugation. 
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Appendix A 
Appendix 


A.l General SU(3) formulae and definitions 


In this Appendix we present the formulae and definitions which will be used 
in the following calculations. 

The Gell-Mann matrices are: 


Ai = 

A 4 = 

A 7 = 


0 1 0 A 

100 

000 / 

0 0 1 A 

000 

100 / 

0 0 0 A 
0 0 -i 
0 * 0 / 



H 0\ 

00 a 3 = 

0 0 / 

0 -i A 

00 a 6 = 

0 0 / 

/ 1 0 0 
0 10 
\ 0 0 -2 


1 0 0 A 

0-10 

000 / 

0 0 0 A 

001 

010 / 


(A.l) 


The actual generators of the 577(3) group/algebra we take as the half of 
the lambda-matrices 

T a = ^A a (A.2) 

They satisfy (for any SU(N) group) the following commuta¬ 
tion/anticommutation relations: 


[T a T b \ = ifabcTc 


(A.3) 


{T a T b } 


J^^ab^N + 4 d abc T c 


(A.4) 
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with f a bc being the absolutely antisymmetric and d a bc - absolutely symmetric 
thensors of the third rank. Traces of the generators and their products satisfy 
the following relations: 



TrT a = 0 

(A.5) 

TrT a T b = i 5 ab 

(A. 6 ) 

TrT a T b T c = | ( d abc + ifabc ) 

(A.7) 

TrT a T b T a T c = ^ N $bc 

(A. 8 ) 

Also 



-j^ 

TaTb = 2 

| 

~J^d ab T (d abc T ifabc)T c 

(A.9) 

r pij r pkl _ 
a a ~ 

= 2 ~ 

(A. 10) 

Structure constants of the group satisfy the Jacobi identities 

fabc feed febefaed fdbeface 0 

(A. 11 ) 

fabc feed febefaed H - fdbeface 0 

(A. 12) 


From the antisymmetric nature of the f a bc directly follows 

fabb = 0 

for any b. Similar self-contraction relation holds for d a b c , where we put the 
explicit summation to distinguish it from the previous formula. 

^ ^ dabb 0 
b 

The contraction relations for / and d are 


faedfbed N8 a b 

(A.13) 

N 2 - 4 

(A.14) 

dacd^bcd ^ab 
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A.2 One Link Integral 

We will use here the Gell-Mann parametrization of the arbitrary SU(3) matrix: 


. OLA X A 

U = e 




e 2 


(A.15) 
(A.16) 


where A* are the Gell-Mann matrices (|A.1|) . For small cp we can expand 


.(Xi\i OLj i OLj\i\j ,OLj i OLjOL} i \j i \j\} z 

U ~ 1 + i ---- — i ---— 

2 8 8-3! 

x i OLj i \i OLiOLjXiXj .OiiOLjOLkXiXj Xk 

— -— - + i- 


(A.17) 


8 


8-3! 


Let us calculate the measure first. Following the standard procedure we 
calculate differentials: 

.dotiXi (doiiOLj + oiidoij) XiXj . idoiiOijOik + oiidoijOLk + cpoodo^) XiXjXk 

du ~ i -----—-- — ------------— 

2 8 8-3! 

.ciccjA* doiiQtj (AjA j -I- A^Aj) . doiiCtjCt^ (AjAjA^ -I- A^AjA^ -I- X^XjXi) 

“ l ^2 8 1 8~3! ““ 

.daiFi doiiOijFij .daiaja^Fij]^ 

,_ 2 8 * 8~3! 

I .d^Fi da^jFij .daiajaj^Fijj^ 

C = 8 +1 8G3! 

Now we calculate the ’’metric tensor” in this parametrization: 

m m / da.nFi danCLiFi ,• daiOiauFauX 

TrdUdU ] ~ Tr i—— - '_LJL _ % — ' * ljL x 

2 8 8-3! / 


x 


■ da a F a da a a b F ab ,da a a b a c F abc 

-i -b i - ; - 

2 8 8-3! 


da^i da a F a daiajFij da a a b F ab daiFi da a a b a c F abc 

= I r — ---1--— -2- 


2 2 
— 9 iadaida a 


8-3! 


(A.18) 


with 


TrFiF a ajabTrFijFab a b a c TrFiF abc 

9ia = —xx-h - — -2- 


2 2 ' 2 2 (2 x 2!)^ 

x.jui b j. vF ij j 


2 2 x 2 2 • 3! 


dja a j a b T r Fjj F ab a b a c FvFiF ab c 

' T + 64 8 G3! 


(A.19) 
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Qil 2 Cjn 1 -f- OijOLk 


(TrFijFik TrFiFi 


Ikj 


24 


V 32 

According to [Q the invariant measure is then 

d/i = a/ det g ]4 dec* 

i 

Using the following representation of the determinant 

det q a = exp Tr log qu 


(A.20) 


(A.21) 


f TrFijF, k TrFtFitj 

log*, = 0 , 0 *^—- — 


rr , (TrFijFik TrFiF ikj \ _ 

Tr log qu = ^ a.ja k I -—-—- ) = aja k T jk 


where we defined (and calculated) 


rp _ / TrFijFi k TrFiFi k j\ _ 1 

F — / j 1 qo oa I ~ a °F 


(A.22) 
(A.23) 


32 24 y 4 

So the expression for the measure takes the following simple form 
d/i = \J det g ]4 dai = CeF» ai ° l J^dct* 


(A.24) 


We also define 

Gij k i = TrXi\j\ k \i 

Explicit calculation leads to: 


Gij k lOLiOLjOL k OLi 2 ^ ^ 


2 2 
a i a j 


hj 


-'X a i a i) 


Now let us define the following average 

/ OO 

J4 da * A exp 

-OO 

We will need these averages: 


77|-4 

(1) =^; (ay.) = (a?) =% 

w 7 /y 4 5 \ 1 3 / 7 4^/6 ’ \ l /'y 4 ^y 6 


(A.25) 
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7 


(A.26) 


E 


hj 


a^cpj 


= ( 8 * 8 - 8 ) ( 




+ 8(af) = 20 


7r 


^/O 


Now we can calculate the statistical sum: 

Z — f d^A Tr ( u+v ''> = e» f d^-A( T ^-v)+ T ^-u<)) 


P ( a i OL j a k a l'^ r ^i^j^k^l 

8Af! 


_ r _ /3 /«*«* 

Cer / e ajakTjk e 2Ar v 2 

„ /* /_ r/3r ,1x1 P a i 0 t j 0 i k a l G ijkl\ 

Ce^ / e (^U^ + S^J 2N 8-4! J _ 


Defining the new coupling 


P 1 

7 “ * AN + 8 


and expanding the integrand we get 
Z = Ce p [ e~ aiaaSi - 


/3 Gijki 

1 + 7TT7 - - 7T a i a j a k a l 


2N 8 • 4! 

Now we will use the formulae for the averages (|A.26|) and get 

A P 20 


Z = Ce?- 


71 

(J- + I) 

V AN ~ 8 / 


1 + 


2iV 96(^ + |) 2 J 


(A.27) 


(A.28) 


(A.29) 


(A.30) 


At this point we set N = 3 (in fact our derivation is not for arbitrary N at 
this point, but we kept it variable for transparency) and obtain for the string 
tension 


a — — In 


= — In 


d In ( Ce p n 4 

(4+1) L 


1 I £ 20 

e 96(4+irj 


dp 


d (j3 — 4 In (2 p + 3) + In ^1 + 20^^^ j 


dp 


- - In 1 - 


2/3 + 3 


+ 


20 _ on £ 

( 2 / 3 + 3) 2 ( 2 / 3 + 3) 3 


1 + 20 


j9 


(A.31) 


( 2 / 3+3 r 


and expanding to the second order in 1//3 we arrive to the final answer: 

<t = | + 4 + 0( 4 ) (A ' 32) 
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Figure A.l: Feynman diagramms for the second order term in the potential 


A.3 One-loop Feynman diagrams for the vac¬ 
uum polarization 

Here we calculate the Feynman diagrams for the n^Q. We define the following 
variables: 

Q • % 
q^ = 2 sm — 


2 

% 


= cos 2 


-2 -2 , -2 - - 
q =q + % = 2 ^ %% 


m =0 


? = 




m= 1 


Feynman rules we need are 


(A.33) 
(A.34) 

(A.35) 

(A.36) 


1. 3-line vertex: (a,p, p) ( 6 , q, u) (c, k, A) gives factor 

igf abc - q)^c Pv + 5 Xp C P - k) v c qx + 8 vim (q - p) x c kfi 


(A.37) 


2 . propagator (a,i)[q](b,j) gives factor 






(1 — <5g 0 ,o) + a 



(?) 



q 2 


(A.38) 
(A.39) 
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Non-static part is 




(A.40) 


EC (?) = 0 


(A.41) 


First diagram 

The first diagram (Fig. |A.l| , left) has two internal lines and two vertices, 
that is 

(a, q, /i) (d, q' - k, A) (c, -q', p) and (c, q', a) (d, q - q', r) (6, -q, u) 

Hence, according to the Feynman rules above we have 


Q i (?) = 


L 2 L( 


E 




D% (?' ~ q) D* (q') x 


I C n > 
' P 


x f adc [<Ap (-?' + ? - ?% c 9A + 5 W (q + q') x c {q ,_ q)p + h A p (?' - 2g) 

X f cdb [5 rv (-q + q' - q) a c g ; + (q + q') T c {q r_ q) ^ + 5 ar (q - 2 q’) u c qa \ 

(A.42) 


x 


Now we set /i = z/ = 0, and we calculate for q = 0 

«-<°> = zfe£(vfe) 

x f adc x (-2g') 0 + h 0 p (?0 a c 9o + ^ A0 (^)p c 9 o 
x f cdb [h T o (?% c q' 0 + ^o-o (?0r c q' 0 + ( — 2g') 0 ] (A.43) 


Only the first term in the first bracket and the last term in the second survive 

Qi (0) = E ( y=) (,') £>“ (,') (-2 9 ') 0 (—2 9 ') 0 (A.44) 


Now we insert the propagators and omit the prime 


Q i (o) 


A?L, 


E 



fadejedb / ^ 





q 2 


QjQj 1 \ 2 
q 2 ql) 


- sin 2 g 0 
(A.45) 
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Using identity (|A.13|) and setting a, b = 0 


N 

Qi (0) = c 


D—l 


EE 


L 2 L 0 4-; V J2c~ f J \q 
k,k 0 ^ov =1 ' v 1 ' 


4(4 


QiQj\ , QiQj 1 \ 1 • 2 

7 sin Qo 

q / q QoJ 4 

(A.46) 


Ah 


D 1 {mj + Mo) 21 ■ 2 
sm g 0 


2L ' L0 + % 2 ) 2 Qo 4 


Ah 


\ ^ 1 . 2 \^ QiQj 4" ^QiQj^ijQo + Ajt/o 

Z, 7 sin r 4o 2^ - 


2L 2 L 0 ^ 4 

9,90^0 ij =i 

4c2 _ at, 


Ah 


2L 2 L, 


E? + 2 ttf:E 


(<f + g 2 ) 2 $ 
Mo 


= - 2 - 


Ah 


<?;Q07^0 


E 


« i2i °« M 2 + © 


cSs„ 2 Av v- 2 A-.!I„- li 


EEh+ 


2Ln 


(A.47) 


L * L ° 7<2o^o (^ 2 + Qo) L ° go^O 

Second diagram 

Now let us calculate the second diagram on (Fig.|A.l|, right). For Feynman 


rules see 46 


<w°) = Ae (A)b“(*/) 2 e/™«/ o 


L 2 L 


° k' o ^0,k’ ^ 2Cf 

N r . 


^Or^OpCk'pCkg A &p T C2k' 0 g^OrAAp 


— ^OpKK + i2^0r<4op M) 


L 2 L, 


o J 


E 


k’ o ^0,k' 


N c 

L 2 Ln 




(A.48) 


/CQ7^0,fc 


2 ko . 2 ^0 0 • 2 ^0 1 > 2 

C 2 kn = COS K 0 = COS-Sill — = 1 — 2 Sill — = 1 - k n 

Zk ° u 2 2 2 2 0 


(A.49) 


N 

Qa( 0) = iVc 


l JV C 1 _ po At c 1 

L 2 Lp ^ 2 + L 2 L n ^ t-2 i 72 + l2t 2-*t I 2 

k,ko^O k,ko^O 0 A:,/co7^0 0 


A4 \ ^ 1 |fcp A r e 1_ 1A4, _ i\ 

L 2 Lq ^ t- 2 _i_ £.2 Eg k 2 2 Lq 

u fc,fe 0 ^0 ^ w «, 0 u fc 0 ^o u 


(A.50) 
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In addition to these two diagrams we need the corresponding ghost diagrams 
and the diagram for the measure. All of such were calculated as in |29| and 
amount to the term 

N " ^ 1 (A.51) 


L ° ko^o K 


Adding all these we get fl 


ns 

00 


rioo(o) = 




EE 

s fco^O k 


9r 2 k 2 
zc fc 0 K - 0 

'k 2 +B 


1 - fcg/2 

k 2 + k 2 


(A.52) 


nSo(o) = 


LnL 2 


EE 


2 cos 2 f4sin 2 f 


2 sin 2 f 


k 0 ^o kl ,k 2 K ( cos2 + cos2 k2 + cos2 

The infinite volume limit is obtained by replacement 

'2tt\ 2 


L s ) 


E 


d z k 


2 cos 2 k\ + cos 2 k 2 + cos 2 k 0 
(A.53) 

(A.54) 


ftoo( 0 ) = j- E 

We dehne the following 


r d 2 k J 

f 8 cos 2 f sin 2 f 

1 - 2 sin 2 f | 

L 4vr 2 | 

1 (cos 2 k\ + cos 2 fc 2 + cos 2 k 0 ) 2 

cos 2 fci + cos 2 k 2 + cos 2 k 0 f 
(A.55) 


, • 2 A™o\ 

( = 8,11 UrJ 

o( f \_ 1 r _ dk\dk 2 _ 

47T 2 E 2 (1 + f) — cos (fcx) - cos (fc 2 ) 

and IIqq can be presented in the more convenient for the practical calculation 
form 

n«°> = i E H (1 - a) SW -»(1 -0 ■(^) 

77-0 — 1 

(A.56) 
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In a more compact way it can be expressed as 


w = nss(0) = E Yx G( x) ' x = 1 ± ' (a ' 57) 

0 no=l ' 0 

G(x) = sinx cosx5'(sin 2 x) (A.58) 

In the integral defining S we replace the inverse denominator in a following 
way 


S(t) — / dkidk 2 / dyexp[—y (2 (1 + t) — cos ki — cos k 2 )\ (A.59) 

J-TT Jo 

and identify two modihed Bessel functions Iq to get 

poo 

S(t)= dyexp[-y(2(l + t))]I$(y). (A.60) 


This can be expressed via hyper-geometric function |42 


as 


S(t) = 


HO 


1/2 


-F (1/2; 1/2; 1;*), * = 


t (t + 2) 


This ansatz relates the behaviour of S (t) at low t to that of F at large negative 
z. Using (421 again one finds 




n =0 


Tin+ 1/2) 
T (n + 1) 


z~ n [log (-z) + 2+(n+l)~ 2+ (■n + 1/2)] 


(A.61) 

A first approximation Wo to lUis obtained from the standard relationship be¬ 
tween a sum and an integral: 


Wq = — 

2i r 


/ dG(x) + - 

'■k/Lo ^ 


'dG(x) 

, dG(x) 


dx 

X =T~ 

L 0 

dx 

X=7T—-JJ- 
L 0 J 


= + 1 dG ^ 
7T L 0 2 dx 


(A.62) 
(A.63) 


' L 0 


where use has been made of the symmetry of dG J^ under x n — x. Up to 
terms which vanish as 7t/Lq —» 0, Eq. (|A.61|) gives: 


Wq ~ 


2vr 


3 log 2 


+ 1 - log 7T + log L 0 


74 





























For that part of G(x) which is analytic at x — 0, this is the right answer. 
A correction to the constant term in L 0 however arises from the logarithmic 
singularity of G(x) (Euler-Mac-Laurinformula, P2f ). As x —> 0, we have 'jjj ~ 
— i log up to analytic terms, whose contribution W l ° 9 to W from Eq. ( |A.61| ) 
can be computed using 


Lo/2 

X] log ^ = log 

—1 -^0 

n 0 =l 


(io/2)| 

T Lq/2 


~ log 


(2e) L °/ 2 


and found to be 

W l ° 9 = ~ [log(2L 0 vr) - L 0 (l + log2)]. 

Ztt 

This we compare to the contribution Wq° 9 to Wq of the same singular part of 
namely 

K 9 = ~ [log(2L 0 ) - L 0 (l + log2)], 
so that in the limit Lq —s► oo the final result is 


IE 


Wo - Wn° 9 + 


'o 

1 

2n 


- log 2 - 1 + log L 0 


(A.64) 
(A.65) 
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A.4 Free field theory 

We present here a solution of the free-held theory, used in Chapter 2 for check¬ 
ing simulation validity for the Higgs sector of the 2D model. Statistical sum is 


Zn = / e s ° 


(A.66) 


with action being 

P 


h 


Tr A 0 (x)A 0 (x + i) — 2 Tr Al(x) — — Tr A 2 0 (x) ) (A.67) 


X \ 2=1 


x= (i,j) i = 0,...,n x -l, j = 0,... ,n y - 1, 


(A.68) 


Fields belong to the 577(3) algebra and can be expressed via generators 


Aq(x) 71“ (x) real 


Using the tracelessness of the generators, we get 

Tr Ao{x)A 0 (y) = -]-^A a 0 (x)A a 0 (y) 


We may rewrite the action in the following form 

2 

^ A^lnA O 4 

■ Hi 


h 


a 'x ) 2 - —An(x ) 2 


s » = ~ E E E +*) - 2 Av, 2 

a x \ 2—1 

=-ivEE(E {Aq(x) — Aq(x + i )) 2 + hAo(x 


a x \ 2=1 

and diagonalize it using the Fourier transform: 

1 

fly £ 


^) = — 


ik-x 


/ 2tT 27T \ / \ r\ / \ 

k = (— q, — r) g(r) = 0,. .., n x (n y ) 


n-r n„ 


In Fourier space 


E* 


x) 2 = 


^x^y 


EwU 


(A.69) 


(A.70) 


(A.71) 
(A.72) 


(A.73) 
(A. 74) 
(A.75) 
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J2(A(x + i)-A(x)) 2 =^J2E A(k)(e ik '~ l - 1 )e ik - x A(k')(e ik ' i - l)e ik ' x 

x n x n y x kk : 

J2 e i {k+kY ' X = n x n y 5 k ,- k ' (A.76) 

X 

A{-k) — A(k), —k = 2n — k (A.77) 

£ (A(X + 0 - ^M) 2 = -V £ \m\\e ,k4 - l)(e- ifc ' S - 1) 

x n 'x n y k 

= — Y\A(k)\ 2 (l - cos(k ■ i)) (A.78) 

n x n y ^ 



77 



d log Z 0 
dh 


--(N 2 - 1) V--- 

^ 2^(1 —cos (k-i)) + h 

n x — l n y~l 

{N 2 - 1) V V _-_ 

9=0 r=o 2 (l - cos(gg)) + 2 (l - cos(gr)) + h 


/ 7vr9 i \ at- n x 1 n y -• 

(Tr v4 2 (a;)) = ^ ~ V V _I_ 

^ ^ 2 (l - cos(gg)) + 2 (l - cos(gr)) + h 

(A.84) 

This quantity can be computed numerically for any given lattice and is very 


suitable to test the simulation of the Higgs sector of the model. 
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A.5 Weak field expansion for Polyakov loops 


One may also test the measurement of the Polyakov loops correlations using 
the small held expansion 

Po,x = ( Tr exp ( A 0 (0) Lq) Tr exp (A 0 ( x ) L 0 ) ) (A.85) 

Po, x = I n dATr exp (A (0) L 0 ) Tr exp (A 0 (x) L 0 ) exp S 0 (A.86) 

X 

We expand Polyakov loops, assuming that A 0 is small: 

exp (A 0 (0) Lq) = I + A 0 (0) L 0 + —Aq (0) Lq + .. (A.87) 

Pq.x = (tt + Aj (0) Lq + —A 0 (0)" L^j Tr + A 0 (. x) Lq + -A 0 ( x ) 2 L^j ^ 

(A.88) 


= <!>/) + l -Ll (Tr (4, (0) 2 )) + \l\ (Tr (4„ (x))' 2 ) 
+ ht(Tr(A 0 (0f)Tr(A o (xf)) 


(A.89) 
(A.90) 


TrA o (0)A 0 (x) = —-^^Ag(0)Ag( 


(A.91) 


F(0, x) = (Tr (A 0 (0) 2 ) Tr (A 0 ( x ) 2 ) > = -- A“(0) 2 £ A b 0 (x) 2 j (A.92) 

To calculate such correlation function we proceed in a standard manner, 
introducing the sources in the action 

S[J| = 41 ArT. -E 4 ^Renm(k) 


k \ a,b 


N n x n y 


A ab(k) = ((1 - cos (A; • ?)) + h) 1 5 a ,b 
New partition function is 

Z[J\ = Z 0 exp 

V x y k a,b ) 


(A.93) 
(A.94) 


(A.95) 
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and the required correlation can be calculated taking the derivatives with re¬ 
spect to the sources and then setting them to zero 


F(0, x ) 


F( 0, x) = 


4 \ n T n 


x'°y 


ex P (ip'R) ( A o(p)f {Kip ')) 2 


a,b,p,p' 


* exp 1 E £ nk)A at (k)j\-k) 

4 V E“P<^) 

x y/ c,d,p,p’ 


X 


d 2 


d 2 


j c ( p'Y J d {p) 


1 (3 1 

2 eX P I o 


2 N n x n y 

* k a,b 


]T]Tj a (fc)A ab (£;)J fe (-<fc) 


l 2 

4 N 2 ( n x n y ) 4 


exp (ip R) 

c,d,p,p' 


(A cc (p')) (A dd(p)) + <W (2A cd (p')) 2 

(A.96) 
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A.6 Factorization of Matterfields 

First, we have to derive some identities which will be used to prove the factor¬ 
ization of the higher-order averages. Let A be an arbitrary matrix belonging 
to the SU(3) algebra. Such matrix can always be transformed to the diagonal 
form. Then on diagonal will be its eigenvalues, which we will call 71, 72 and 
73. Tracelessness implies 73 = —71— 72- Our consideration is based on the 
general expression for the determinant 

\ogDet(l — tA ) = Tr log(l — tA) (A.97) 

Let us calculate the left-hand side of it: 

Det{ 1 - tA) = (1 - £71) (1 - £72) (l + t (71 + 72)) 

= 1 ~t 2 ( 7 1 + 72 + 7172 ) + t 3 ( 7 x 72 + 7i7 2 2 ) 

The quantities which we are looking for are: 

A 2 = TrA 2 = 7 x + 72 + (71 + 72) 2 
= 2 (Ti + 7 2 2 + 7i72) 

A 3 = TrA 3 = 7i + 72 - (71 + I 2) 3 = 

= -3 (7?72 + 7172) 

and we have 

Det( 1 — tA) = 1 1 2 -^ t 3 

Now we shall expand logarithms on the both sides of the identity. As long 
as this expansion should be valid for arbitrary values of t, we can match the 
coefficients near any power of it. Left-hand side: 

\xiDet(l-tA) = --A 2 t 2 --A 3 t 3 --A 2 2 t 4 --A 2 A 3 t 5 -—Alt 6 -—A 3 2 t 6 + 0 (t 7 ) 
v ) 2 2 3 3 g 2 6 2 3 lg 3 24 2 t- v j 

Right-hand side: 

Tr log(l—f A) = --t 2 TrA 2 --t 3 TrA 3 --t 4 TrA i --t 5 TrA 5 --t 6 TrA 6 + 0(t 7 ) 
’ 2 3 4 5 6 v ’ 

and we immediately get: 

A 2 = 2 TrA 4 (A.98) 

-A 2 A 3 = TrA 5 (A.99) 

6 
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(A.100) 


-Al + = TrA 6 

3 3 4 2 

Mean Field Calculation for A 2 2 : 

Here we will work within the mean field approximation: 

A a (x)Ap(x) -> ^5 ag A 2 (x) 

2 A a {x) = —^2A a (x)A a (x)Ap(x)Ap(x) 

Using Wick’s theorem we get: 

(2 A a (x)) = ( A a(x)A a (x)Ap(x)A p (x)) = ^ ( A <*(x)A a (x)) (A p (x)A p (x)) 

a,f3 a, (3 


+ 4 2 ( A a(x)Ap(x)) (A a (x)A p (x)) 

a,0 



{A n (x) A m (0)) = A nm (x) + (An (x)) (A m (x)) 
And we arrive to the final result 

( A \ (x)) = A 2j 2(0) + A^ 

Calculation of A 33 

( A 3 ( x ) ^3 (0)) = A 33 (x) + (A 3 (x)) (A 3 (x)) 


In our phase 


(A 3 (x)> = 0 


(A.101) 


so 

A33 (0) = (A 3 (x) 2 ) 

TrA 3 = ^2 AiA.jAkTrTiT.jTk = ^ AiAjA k ^ ( d ijk + if ijk ) = ^ ^ dijkAiAjAk 

ij k ij k ij k 

where we used the antisymmetric properties of 

(TrA TrA ) -- EE dijkdlmnAiAjA]^AiA m A n 

ij k Imn 
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Now we will use the property of d - it being zero after contraction when two 
indices are the same and its symmetry which leaves us only with six equal 
terms 


11 (l \ 3 IV 2 — 4 


(TrA 3 ) = jj(jA 2 ) * d ijkdijk = 77 ( 7^2 I 6 


ijk 


44 V 4 


N 


(■ N 2 - 4 = 


Or 


(TrA>) = ^ 
^ 33 (°) = ^ A 2 


Calculation of (A 2 A 3 ) 


11 


t^HiiEE A n A n AiAjA k d. 


ijk 


n ijk 


11 


5i 3 EE ( A n A n Ai ) (AjAk) dijk+ 


24 
11 


n ijk 


#EE {A n AiAj) ( A n A k ) dijk+ 


24 

11 
24 


n ijk 


EE (blnbln) ( AiAjAk) dij k 


n ijk 


(A 2 Az) — 

4 _ _ ^ 

= 24 ^ ^ ^ ^ 1 ( A n A n Ai) -A 2 8jkdijk+ 

n ijk 

+ M 6 EE (A n AiAj ) ^ A 2 5 n kdijk~^~ 

n ijk 

+ Mee (AiAjA k ) dij k 

n ijk 

The first term falls out completely - contraction of d with two equal indices. 

I I _ _ I 

(A 2 A 3 ) == 24 dijk^nk (-A-nAiAj') ^A 2 ~\~ 

n ijk 

+ 5iEE ^dijkA 2 (AiAjA k ) 

n ijk 
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Summing over n : 


{A 2 A 3 ) — - -6-A 2 ^ d ijk (A k AiAj) + 2 4^4^ 2 ^ij k ( AiAjA k ) — 

ij k ij k 

= 24^ + d'A (AkAiAj) = --(6 + 8)A 2 A 3 = -A 2 A 3 

ijk 
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Correlation Functions 


Now we turn to the correlation functions. Square brackets denote contrac¬ 
tion and we will silently use the formulae above. 

(An (x) A m (0)) = A nm (x) + (A n (x)) (A m (x)) 


A 24 (R) = (A 2 (. R ); kl 4 (0)) = - <4L 2 (R ); A* (0)> = 


1 
2 

N 2 + 1 


TV 2 4-1 

4 (R) > (°) 


= (*) ^2 (0)]] [[A (o)]] 


5 


N 2 _ l ^ 22 ( R ) A 2 - -^A 2 2 (. R ) A 2 


A 35 (. R ) — (kL 3 (i?); 4L 5 (0)) — - (A 3 (i?); A 3 (0) A 2 (0)) — 

o 

= g[[«A (*)>>; «A(o) A (<>)>>]] = 

A AT 2 4- c; 

= [[<<^3 (R))) ; «^3 (0))) i(A 2 (0)))]] = 

A /V 2 4- 5 

= g jy 2 _ i (t[^ 3 (R) ; A 3 (0)]] [[A 2 (0)]] + [[A 3 (. R ); A 2 (0)]] [[A 3 (0)]]) = 
= -a 33 a 2 

24 

3 44 (fl) = (A (J?); A ( 0 )> = [[((A («)»1 <(4 (o)})]] = 

= G^t) 2 [W W4l(°)]] = 

= tlAl+it 2 1 2 [[A (R) ; A (0)]] [[A (R) ; A (0)]] + \ 
\2N 2 -l) \ +4 [[A 2 (/?); j4 2 (0)]] [[A 2 (0)]] 2 J 

= (2^22 ( R ) 2 + 441 22 (R) Al) = 

= (2A 22 (R) 2 + 4:A 22 (R) Al) 
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(R) = Ps (fl) ;A S (0)> = [[«A 5 (if))) ; (P 5 (0)))]] = 

= (g) (¥r=i) p2 (R) A, (R ); A 2 (0) A, (0)]] = 

[[A 2 ( R ); A 2 (0)]] [[A 3 (R) ; A 3 (0)]] + \ 

+ [[M (R ); ^3 (0)]] 2 + 

+ [[A 2 (R ); A 2 (0)]] [[A 3 (-R)]] 2 + 
+2[[A 2 (R)-A 3 m][[A 3 (R)]][[A 2 (R)]\ + 

+ [[^3 (R) ; A 3 (o)]] [[A 2 (i ?)]] 2 / 



/ 5 \ 2 / N 2 + 5 \ 2 
V6y Viv 2 -iy 


In the last step we used the symmetry relation (|A.101|) 
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